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ABSTRACT 


The aeroelastic flutter a nd d ivergence behavior of 
rectangular, graphite/epoxy/cantilevered plates with vary¬ 
ing amounts of bending-torsion^stiffness coupling is inves¬ 
tigated for incompressible flow. A general Rayleigh-Ritz 
formulation is used to calculate flexibility influence coef¬ 
ficients, static deflections, divergence velocities, vibra¬ 
tion frequencies, and flutter velocities. Flutter calcula¬ 
tions are done using the U-g method. Test plates were con¬ 
structed and subjected to static, vibration and wind tunnel 
tests. Wind tunnel tests indicated static deflections, 
divergence instabilities, bending-torsion flutter at low 
angles of attack, and stall flutter at high angles of attack. 
Bending stiffness and first bending frequencies showed good 
agreement between theory and experiment. Torsional stiffness 
and first torsion frequencies were not accurately predicted. 

”by the theory for highly coupled plates. Divergence velocities 
and reduced flutter velocities showed reasonable agreement 
between theory and experiment. Test plates with varying 
amounts of coupling exhibited markedly different stall flutter 
characteristics. „ 
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CHAPTER I 
INTRODUCTION 


1.1 Background 

Aeroelastic tailoring is broadly defined as the tech¬ 
nology to design a lifting surface which exhibits a desired 
aeroelastic response. Desired aeroelastic responses most 
often considered are maximization of flutter and divergence 
speeds. Other aeroelastic responses can be important de¬ 
pending on the application. They include control reversal 
speed, camber changes as a function of load and speed, and 
angle of attack changes as a function of load. The aniso¬ 
tropy of advanced composite materials or, more specifically, 
the designers' ability to control that anisotropy by selec¬ 
tive lamination, makes it an attractive material for aero¬ 
elastic tailoring. 

Aeroelastic tailoring, which exploits the advantage 

of advanced composites, has received considerable attention 

in recent literature. N.J. Krone, Jr. 1- concluded that 

forward swept wings without divergence or weight penalties 

may be possible through the use of selectively laminated 

2 3 

advanced composites. T.A. Weisshaar ' extended, analyti¬ 
cally, Krone's conclusion to potentially practical wing 
designs. Weisshaar concluded that the binding-torsion 
stiffness coupling of anisotropic advanced composite 
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materials was the key to eliminating divergence in forward 

4 

swept wings. V.C. Sherrer, T.J. Hertz, and M.H. Shirk 
conducted a series of wind tunnel tests using simple plate¬ 
like models of a forward swept wing. These tests essentially 
verified Weisshaar's conclusion, and also showed that existing 
analytical techniques (computer programs) would adequately 
predict the divergence dynamic pressures for most test con¬ 
ditions . 

Under contract to the Air Force Flight Dynamics Lab¬ 
oratory, General Dynamics Corporation developed a wing 

5 

aeroelastic synthesis computer program (TSO) . The computer 
program was intended to be a preliminary design tool for 
optimization of wings with advanced composite structural box 
skins. To this end, it used a direct Rayleigh-Ritz energy 
formulation to model the structural deflections. The prog¬ 
ram was capable of optimizing a wing skin design for several 
different constraints, simultaneously. Published aeroelastic 
tailoring studies ' ' using TSO have not examined, in depth, 
the effect of structural bending-torsion coupling on flutter 
speed in the absence of changes in other variables. Addi¬ 
tionally, only very limited wind tunnel data are available 
which examine the effect of structural bending torsion 
coupling on flutter speed, or even that verify the accuracy 
of flutter speeds generated by TSO. Finally, the phenomenon 
known as stall flutter has received only limited attention 
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8 9 

as it relates to metal lifting surfaces ' , and virtually 
none as related to advanced composite lifting surfaces. 

The references cited in this section by no means com¬ 
prise a complete list of work done in this area. They should 
be viewed more as an indicator of where emphasis has been 
placed. One obstacle to conducting research in aeroelastic 
tailoring is that much of the work done by the larger aero¬ 
space corporations is proprietary and consequently is never 
reported in literature accessible to the public. 

1.2 Objectives 

This study will attempt to ascertain the effect of 
varying amounts of bending-torsion coupling on both the 
divergence and flutter speeds of an unswept lifting surface 
in incompressible flow. The lifting surface will be ideal¬ 
ized by a cantilevered, rectangular, flat plate constructed 
of laminated graphite/epoxy, and having a half-span aspect 
ratio of 4 (8 by standard aerodynamic convention). 

Specific objectives are to 

(1) Develop an analytical formulation to predict static 

deflections, natural vibration frequencies, divergence 
speed, and flutter speed for idealized lifting surfaces. 
The formulation must be applicable to a plate having 
substantial bending-torsion coupling, but limited to 
a mid-plane symmetric lamination arrangement. 













(2) Evaluate the effect of bending-torsion coupling on 
the idealized lifting surface by performing static 


deflection, vibration, and wind tunnel tests on sel¬ 
ected cantilevered plates having varying amounts of 
bending-torsion coupling. 

(3) Investigate the effects of bending-torsion coupling on 
stall flutter velocity. This will be limited to an 
experiment evaluation. 

(4) Determine the accuracy of the analytical formulation 
by comparing theoretical and experimental results. 

1.3 Organization 

Chapter II develops a general Rayleigh-Ritz energy 
formulation to model lateral deflections of a laminated 
cantilever plate. Solutions for both two term and three 
term deflection equations are derived. The two term general 
solution is than applied to a static deflection problem to 
calculate a 2 x 2 matrix of flexibility influence coef¬ 
ficients, and to a divergence problem to calculate the 
divergence velocity of the idealized lifting surface. Both 
two and three term solutions are applied to a free vibration 
problem to determine the lower two or three* vibration fre¬ 
quencies. Finally, the two term solution is applied to a 







classical bending-torsion, potential flow, flutter problem 
using the U-g method. 

Chapter III presents the experimental test apparatus 
and procedures for the static deflection tests, free vibra¬ 
tion tests, and wind tunnel tests (flutter and divergence). 
The test results are then discussed, commenting on possible 
sources of error. 

Chapter IV compares the analytical and theoretical 
results, identifies areas of good and poor agreement, and 
comments on sources of inaccuracies in the theoretical 
analysis, as well as ways to improve it. 

Chapter V summarizes the conclusions reached in the 
previous chapters and makes recommendations for further 
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CHAPTER II 
THEORY 


2.1 Anisotropic Plate Flexural Stiffness 

The flexural modulus (stiffness) components of a lami¬ 
nated advanced composite (graphite/epoxy in this study) plate 
are dependent on both the fiber orientation and stacking 
sequence of the individual laminae (plies>. To simplify 
fabrication of test specimens, only laminated plates with 
a mid-plane symmetric stacking sequence are considered . 

The ply angles (0) follow the sign convention in Fig. 1. 



Figure 1. Ply Angle Sign Convention 
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The in-plane, on-axis lamina modulus components (Q^j) 
were obtained from the orthotropic engineering constants 
for Hercules AS/3501-6 graphite/epoxy, from which the test 
specimens were to be fabricated. These engineering constants 
take on different values depending on whether they are ob¬ 
tained from bending or stretching tests. Engineering con¬ 
stants obtained from each type of test appear in Appendix A, 


and their 

validity 

for vibration problems 

is briefly dis- 

cussed in 

Chapter 

IV (Reference 

10 has an 

in-depth treatment 

of this subject.) 

The terms 

are defined as 11 

Q 11 

- e l /(1 

V LT V TL ) 


(2-la) 

Q 22 

= E t /(1 

" V LT V TL ) 


(2-lb) 

Q 12 

- Q 21 - 

V LT E T /(1 " V 

LT V TL ) 

(2-lc) 

Q 66 

= g lt 



(2-ld) 


where v TL = (E T / E L > V LT 

( 6 ) 

The off-axis lamina modulus components (Q^j ) were 
obtained by first defining a set of invariants 11 


X 1 * 


[Q 11 + Q 22 + 2Q 12 ]/4 


Z 2 = 


[Q 11 + Q 22 - 2Q 12 + 4Q 66 1/8 


R 1 = 


[Qll - Q 22 ]/2 


R 2 = 


(Q 11 + Q 22 - 2Q 12 “ 4Q 66 1/8 


(2-2a) 
(2-2b) 


(2-2c) 


(2-2d) 
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and using 

relations 

the invariants in the 

11. 

following 

transformation 

Q u 

s 

I 1 + I 2 + R 1 cos20 + 

R 2 cos40 

(2-3a) 

q<|’ 

= 

I x + I - R 1 cos29 + 

R 2 cos40 

(2-3b) 

q<«) 

w 12 

= 

r i “ Z 2 ~ R 2 cos49 


(2-3c) 


= 

i 2 - R 2 cos40 


(2-3d) 

< 

= 

R^sin20 + R 2 sin40 


(2-3e) 

q (6) 

^26 

= 

j R 1 sin20 - R 2 sin40 


(2-3f) 

whe re 9 

is 

the ply angle (Fig. 

1) . 



The flexural modulus (D^Jfor an n-ply laminate with 
arbitrary ply angle orientation is obtained from 


D.. = Z Q* ( ? k) [(Zv _ 3 )/3] i;j = lf2i6 (2 _ 4) 

JC = 1 ^ 

where 

0 k = ply angle of the (k) ply 

= distance from the mid-plane to the upper sur¬ 
face of the (k) ^ ply (positive above mid-plane, 
negative below mid-plane) 

z k _^ = distance from the midplane to the lower surface 
of the (k) th ply 

Flexural modulus values (D^, D 66' °16^ ^ or 
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[+30-/0] . [+45-/0] , and [±45/0] laminates appear in 
Appendix A. The reader should observe that, using the 
equations in this section, the flexural moduli of a [+6 2 /°1 3 
laminate will be the same as those for a [-0 2 /O) s laminate, 
with the exception of a negative sign on D lg and D 2g for the 
latter laminate. 

2.2 General Rayleigh-Ritz Formulation 

The direct Rayleigh-Ritz energy method is a relatively 
simple, straightforward approximation for the plate deflec¬ 
tions, as required for the static deflection,free vibration, 
divergence, and flutter analyses in this study. The Ray¬ 
leigh-Ri tz method also has the advantage of showing the 
effect of the individual variables on the solution more 
clearly than other more accurate methods, such as finite 
element analysis. Whether or not the results of the Ray- 
leigh-Ritz analysis developed in this chapter correlate 
sufficiently well with experimental results is discussed in 
Chapter IV. To simplify computation, the "wing" is idealized 
by a rectangular cantilevered flat plate with uniform thick¬ 
ness. Further, to allow the stiffness properties of the 
plate to be more accurately depicted, an aerodynamic fairing 
will not be used over the plate. 

The Rayleigh-Ritz analysis begins by assuming a deflec¬ 
tion shape for the structure. If only lateral deflections 









(w) are allowed, the single deflection equation, written in 
generalized coordiates is: 
n 

w = I Y-i (x,y)q. (t) (2-5) 

i=l 1 1 

where y^(x,y) is the nondimensional deflection or mode shape 
of the (i)^ mode (which must satisfy geometric boundary 
conditions for a cantilevered plate), and q^(t) is the gen- 
eralized displacement of the (i)^ mode. The generalized 
displacement is a function only of time and has units of 
length. For convenience and consistency, the coordinate 
system in Fig. 2 will be used for all problems. 

The deflection equation is simplified by assuming 

(1) the plate is chordwise rigid when undergoing torsional 
deflections, 

(2) the plate does not exhibit chordwise bending for any 
of the problems to be developed, 

(3) a single term in the deflection equation for each of 
the desired deflection modes (first bending and first 
torsion; or first bending, second bending, and first 
torsion) will adequately represent the actual deflec¬ 
tions of the plate. A further requirement here, to 
insure rapid convergence, is that the mode shape for 
each term must accurately depict the deflection for 
that mode. The two term deflection equation can now 


be written as 













i 




Figure 2. Sign Conventions 










w 


<j> 1 (x) 4 > 1 (y)q 1 (t) + * a < x) V y)q 2 (t) 


(2-6) 


where <^(x) and (y) are one dimensional mode shapes for 

'f'Vi 

the (i) mode. If the first term represents the first 
bending mode, then <J> 1 (y) is simply a rigid-body translation 
($^(y) = 1) • If the second term represents the first tor¬ 
sion mode, then, by the first assumption, <f> a (y) is a rigid 
body rotation (<j> (y) * y/c). Incorporating these into Eq. 

2-6 yields 

w = ! J) 1 (x)q 1 (t) + (y/c) <J> a (x) q 2 (t) (2-7) 

We can apply the same analogy to the three term deflection 
equation, which is written as 

w = < j) 1 (x)q 1 (t) + <J> 2 (x)q 2 (t) + (y/c) (x) q 3 (t) (2-8) 

In Fig. 2b, the elastic axis (EA) is, by definition, 
located at the midchord of the plate (y - 0). Lateral def¬ 
lection of the elastic axis is designated w £ and rotation 
of the plate about its elastic axis is designated a. One 
observes that w E and a are not generalized coordinates and, 
in fact, q 2 in Eq. 2-7 and a do not even have the same units. 
The units of a are radians, while q 2 has units of length. 

To write w^ and a in generalized coordinates (for Eq. 2-7), 
the following relations are used: 


i 
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W E = w y 2*0 * ♦ j.Wq^t) (2-9) 

a = || = («t a (x)/c)q 2 (t) (2-10) 


The importance of this transformation will become apparent 
when aerodynamic forces are applied. Finally, from Fig. 2b, 
Fg is a lateral force applied at the elastic axis and Mg is 
a moment about the elastic axis. 

The strain energy (V) for a symmetric anisotropic lami- 

. 12 
nate is 


V = 


i / 
z 0 


c 

2 


c 

'7 


lD ll<“'xx> 


2D,_w, w, 

12 'xx 'yy 


D 22 ( "'yy) 2 + 


4D,,w, w, + 
16 xx 'yy 


4D_,.w, w, 

26 yy xy 


+ 4D 66 (w 'xy ) 1 dydx 


( 2 - 11 ) 


where a comma denotes partial differentiation. 


2.2.1 Two Term Deflection Equation 


w 'xx - »'i1l + I *1*2 w, xy = 


w, yy - 0 


( 2 - 12 ) 


where: ( )' = d/dx and ( )" = d^/dx^. Substituting Eq. 2-12 


into Eq. 2-11 and performing the chordwise integration yields 
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V = 


qjl ^ o /(** l ) 2 dx] + 


-j CD. . x. 0 tv cc x, -j 

+ q,[ f (<{, n ) 2 dx + —££■ / ( 4 ) ’ ) 2 dx] (2-1 


2D, 


l 2 24 


The kinetic energy expression for the plate is 


T = 


\ / 
z 0 


c 

7 2 

/ m w dydx 


( 2-1 


where (') = d/dt and m is the mass/area. 


m = PGE fc p 


( 2 - 1 ] 


where p G£ is the specific gravity of graphite/epoxy and t 

P 


the plate thickness 


w » + (y/c)4> a q 2 


( 2—11 


Performing the chordwise integration, Eq. 2-14 then becomes 


T = q* l X Q /%idx]+ q 2 ( / ^dx] 


( 2 - 1 ' 


The change in external work (<5W e ) can be expressed as 


<5W„ 


c 

J 


f c ! P z Sw dydx 


( 2-11 


where p z is a distributed lateral load and 


< 5 W = *, 59 , + 


( 2-1 
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Then, by substitution 


6W e = Sq^ + 6q 2 Q 2 


where 


Q 


1 


l 

f 

0 



P^dydx 


and 



c 

1 V 

_/ p_ *; <J> dydx 
c J c z c 

2 


( 2 - 20 ) 


The relationship between the work and energy expressions 
is obtained through Lagrange's equation, which is a statement 
of Hamilton's energy principle; the basic premise of the 
Rayleigh-Ritz method. Lagrange's equation is 13 


d_ 

dt 


3T ) + 3V 

3< *i 3q ± 


3T = Q i (2-21) 

3q ± 

(i « 1,2,...N) 


where: 


3V 
3<3 


- 0 


3V 

3^ 


dc i a 

2q [ -il- / («{,••) 2 dx] + q,[2D., /<J> * <(,' dx] 

A 2 0 * ■ L0 0 a 


1 2 * -Jl I2D 16 / «W xl + 2 ^2 / <♦;> 

+ -w j dx! 


cD 


11 


"dx 












From 

with 
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^ = 2q, [5^ / dx] 

aq L 1 ^0 1 

i 

4r ( — ) = q, [ me / <}>? dx] 

dt 3q L 0 


9T r me - ,2 - ■. 

- 2c ?2 [ 24 J dx ] 

3 q 2 0 


d , 3T . _ r me , ,2 , , 

dt a • q 2 12 n f dx5 

3q 2 0 


Eq. 2-21, the general equations of motion for the plate 

a two term assumed displacement equation are 

i i 

q [ me f $ 2 dx] + q. [D..C / (<j>") 2 dx] 

i 0 1 L xx o 


+ q 2 [ 2D lg f dx] = (2-22a) 


V r! / *1 dxl + <11 t 2D 16 / *i*i dxl 

d ^ll ^ 2 

+ ^ir/ <♦;> dx 

4D 1 ? 

+ —~ / (♦:) dx] = q 2 

^ o u 


(2-22b) 


To make these equations more compact, nondimensional expres¬ 
sions for the integrals are defined below. These integral 
expressions may be evaluated exactly, or using any of several 
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numerical integration techniques, if the integrand is too 
complicated. 


X 3 

= 

1 

z 

2, 

/ 

0 


(2-23a) 

X 4 

= 

1 

£ 

z 

s 

0 


(2-2 3b) 

X 5 

= 

1 

I 

£ 

f 

0 

<j>^ dx 

(2-23c) 




z 



X 6 

r 

z 2 

/ 

0 

dx 

(2-23d) 

X 7 

= 

z 3 

£ 

/ 

0 

(0'{) 2 dx 

(2-23e) 

X 8 

r 

z 

£ 

/ 

0 

(4^) 2 dx 

( 2-23f ) 

X 11 

= 

l 3 

£ 

/ 

n 

($;) 2 dx 

( 2 — 2 4 g) 


0 


Incorporating the expressions of Eq. 2-23 into Eqs. 2-22a 
and b yields the final form of the equations of motion. 


<*1 

$2 


D 11 C 


2D 


16 


[mcU 4 ] + q L [-p- I ? ] + q 2 [ I 6 J = Q 1 (2.24a) 

r me* T i , _ , 2D 16 t , 

[ TT X 5 3 + q l [ ~T T ~ X 6 3 


£- 

. r 4D 66 , D 11 C " - . , 

+ q 2 c£ {I 8 48Dg g Jl 7 I ll ) 3 


EQUATIONS OP MOTION 


(2.24b) 


In Eq. 2-24b, the term 
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48FD^“ "11 

represents the warping stiffness, which is inversely pro¬ 
portional to the square of the plate aspect ratio ( Z/c ). 
For the plates in this study (AR = 4), this term's contri¬ 
bution will be small. 


2.2.2 Three Term Deflection Equation 

The derivation of the equations of motion for the three 
term deflection equation is almost identical to the two-term 
deflection equation derivation. Therefore, only the final 
expressions for V, T, 5W_, and Lagrange's equation will be 
presented. As one would expect, there will be three equa¬ 
tions of motion instead of two. 


Deflection equation: 


w = $ 1 (x)q 1 (t) + <j> 2 (x)q 2 (t) + ^<J> a (x) q 3 (t) (2-25) 


Strain energy: 


'11'* " ? D-i -i C £ 7 7 

V = [ 2~ Q f (< * , 1 ) dxlq l + 1 ~T~ / <4>2 } dxlq 2 


Di I c Z o 

/ . It \ , 


2D 66 * 


+ 1 -IT" ! ( <> Mx + —ir ! U’) 2 dx]q^ 

0 a 0 a 

z z 

+ [ D n d / dx]q 1 q 2 + [ 2 D lg / ^^dxlq-jq-j 


+ ^ 2D 16 ^2^a dx ^ q 2 q 3 


(2-26) 
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Kinetic energy: 


t = qf [ %£ i ^:dx] + q ; [ ^ ; 4>;dx] 


me 


+ I U / 4>‘dx] + q 1 q 2 [ me / ♦ 1 * 2 dx] 


(2-27) 


Change in external work: 


5W e = Q 1 6q 1 + Q 2 5q 2 + Q 3 <Sq 3 


(2-28) 


where: 


i - 

Q 1 ■ o f -£ /2 d * dx 


°2 * 


f J P 2 4>2 d Y dx 


Q 3 = 


0 


f c /2 P 2 &cx dydx 
~2 


Before applying Lagrange's equation, we eliminate terms 

with i i 

f 1 4> _d.x and / 4>" dx 

0 1 Z 0 1 

by requiring that <t^(x) and 4> 2 (x) orthogonal functions. 
Now, using the same nondimensional expressions for the inte¬ 
grals as before, the equations of motion obtained from 
Lagrange's equations are: 
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(2-29a) 


q 2 [ mcl I 12 ] + q 2 f 
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(2-29b) 
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EQUATIONS OF MOTION 


where 


= Jl 2 / <|>" 4 >' dx 

2 r a 


2 

u \ ^ 


i 1Q = z 3 ; ( 0 n ) dx 


12 


1 2 
2 / <t > 2 dx 

X, o ^ 


The functions used in this study for the mode shapes 
<j> 2 , <j> a ) are combinations of sin, cos, sinh, and cosh. 
These functions are known as ''beam'' modes because they 
very closely model the deflections of an isotropic Euler 
beam. The first bending mode shape (<^) and the first 
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torsion mode shape (4> a ) were obtained from Ref. 14 and the 
second bending mode was obtained from Ref. 15. These mode 
shapes, along with numerical values for the nondimensional 
integral expressions, appear in Appendix B. 


2.3 Static Deflection Problem 


The static deflection problem is formulated as an 
analytical model of the experimental deflection tests 
described in Chapter III. For this problem, the canti¬ 
levered plate is subjected to a concentrated unit load and 
moment, individually applied at the 3/4 length point on 
the elastic axis. Since the result of the experimental 
tests was a 2 x 2 flexibility influence coefficient matrix, 
the same result is desired here, to allow direct comparison. 
The 2x2 matrix requirement dictates that we use the equa¬ 
tions of motion derived from the two term deflection equa¬ 
tion. Observing that acceleration (q) for static deflec¬ 
tions is zero, Eqs. 2-25a and b are repeated here with q 
terms eliminated. 


^ 11 C 

^1 ^ ~ ip ^7 ^ + ^2 ^ 


2D 


16 

T~ 


1 - Qi 


(2-30a) 
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(2-30b) 
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The generalized forces (Q^) obtained from Eq. 2-20 are 


Q 1 = W I Z ] 


Q 2 = Va ( )/C 


(2-31a) 
(2-31b) 


where F fc is a concentrated test load and M is a concentrated 
test moment such that 

£ - Z - 

F = / f 2 p dydx and M - / f 2 y p dydx 

T 0 2 T 0 Z 


The test forces and moments are applied at the same loca¬ 
tions and with the same sign conventions as F £ and in 

Fig. 2b. By the relations of Eqs. 2-9 and 2-10, the def¬ 
lections w_ and a taken at x = %Z are 
E 4 

w £ = ) q x (2-32) 

ct - [ 4> a ( )/c ] q 2 (2-33) 


Assuming unit forces and moments, instead of unit 
displacements, dictates that one will obtain flexibility 
inflence coefficients (c^) by solving the equations of 
motion. The matrix of influence coefficients must then be 
inverted to obtain the stiffness influence coefficients 
(k„) . This is generally much easier to accomplish experi¬ 
mentally than to try and obtain the stiffness influence 
coefficients directly, by applying unit displacements to 
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the structure. To aid in conceptualizing the problem, the 
equations of motion in terms of flexibility influences are 
presented. 




* C 11 F T + C 12 M T 


(2-34a) 


a * c 21 F T + C 22 M T 


(2-34b) 


The flexibility influence coefficients are defined as follows. 


°11 E V f t 


M = 0 
T 


C 12 H W E /M T 


F = 0 
T 


(2-35) 


C21= - “"'Iv 0 


c 22 = a / M tp 


V 0 


To obtain the four flexibility influence coefficients, 
Eqs. 2-30a and b must be solved simultaneously twice, each 
time with a different combination of and Q 2 * The first 
combination (F T = 1 and = 0) will be called "Test 1", 
and the second combination (F T = 0 and M T = 1) will be called 
"Test 2". The values for q^ and q 2 obtained from Test 1, 
when substituted into Eqs. 2-32 and 2-33, yield c^ and C 2 ^ 
as defined by Eq. 2-35. Similarly, q^ and q 2 from Test 2 
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yield c^ 2 and c 22 . The final step, if desired, is to 
arrange the flexibility influence coefficients in a 2 x 2 
matrix and invert it to obtain the stiffness influence co¬ 
efficients . 


2.4 Aeroelastic Divergence Problem 


Aeroelastic divergence is a static deflection problem 
such that the aerodynamic forces and moments applied to the 
cantilevered plate are a function of the torsional deflec¬ 
tion of the plate. Divergence is defined as the point at 
which the restoring forces generated by the structure can 
no longer counteract the aerodynamic forces. At that point, 
by the simple linear theory used in this study, the deflec¬ 
tions increase without bound resulting in structural failure 
Actually, as discussed in Chapter III, both structural and 
aerodyanmic nonlinearities come into play at large deflec¬ 
tions and limit the maximum deflections to finite values. 
However, structural failure may still occur. These non- 
linearities, while important, will not be considered in this 
development. 

The static forms (q = 0) of Eqs. 2-25a and b are re¬ 
peated here. 





+ ^7 



Q, (2-36a) 
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(2-36b) 


The generalized forces (Q^) obtained from Eq. 2-20 are 


/ 4> 1 L e dx 


(2-37a) 


l 4 > 


' IT “e 


dx 


{2—3 7b) 


where 


L 


E 


“e 



p z dy 
Y P z dy 


L_ is the lift per unit length acting at the elastic axis 
of the plate, in the same direction as F E (Fig. 2b). 
is the aerodynamic moment per unit length, again acting 
at the elastic axis. 

For this simple static analysis, two-dimensional 
aerodynamic strip theory was deemed to provide sufficient 
accuracy. The lift per unit length and moment per unit 
length are 

L e = qca Q (a 0 + a e ) (2-38a) 


Mj. = qcea Q (a 0 + a & ) 


(2-38b) 
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where q is the dynamic pressure, e is the distance between 
the elastic axis and the aerodynamic center (quarterchord) 
of the wing, a Q is the two dimensional lift curve slope, 
otg is the rigid plate angle of attack and a g is the flexible 
twist of the plate, such that a n + a. = a (Fig. 2b). 

The lift curve slope of a flat plate of infinite 
length (span) is the well known value 2 it . This two dimen¬ 
sional lift curve slope can be empirically corrected to 
account for the affects of finite span by applying the 
correction suggested in Ref. 13. 


3C, 


3a 


. a [ 

0 l AR + 2 


(2-39) 


where 3C^/3a is the corrected lift curve slope and AR is 
the aspect ratio (2£/c) of an equivalent wing, as defined 
for an aircraft wing composed of two cantilevered plates. 

Equations 2-38a and b must be transformed to generalized 
coordinates before they are incorporated in the equations 
of motion. As before, Eq. 2-10 is used to transform a Q 
only. The generalized forces finally take the form: 


Q x = qc2. 
Q 2 = qei 


5 

3a 

3a 


[ a 0 r i + c I 3 1 


[ a Q I 2 + 


c X 5 5 


(2-40a) 


(2-40b) 


where the nondimensional integral expressions are 
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<f 1 dx 

♦c. dx 
*i+c dx 
♦a dx 


Substituting Eqs. 2-40a and b into Eqs. 2-36a and b yields 
the complete equations of motion. These equations are 
presented here in matrix form to make the solution technique 
more apparent. 


D 11 G 


2D 


l 

2D 


— J 7 

3 ' 




16 


I* - & 


3C, 


3a 


4D 


D..C 2 3C„ 

— (Ig + ——- I,,) 1 - q !* — I- 

cZ 8 4 8D 66 5, 2 11 C 3a 5 




l * ) 


‘ 3C 0 ' 
qcJl 




q el 


3a 

ffl 

3a 


) a r 


(2-41) 


Classical divergence is normally evaluated at zero 


initial angle of attack (a^ = 0), which is also assumed 
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here. Therefore, a unique solution to Eq. 2-41, with 
aQ = 0, can be found by setting the determinant of the coef¬ 
ficient matrix equal to zero and solving for the dynamic 
pressure (q). The dynamic pressure which makes the deter¬ 
minant go to zero is the divergence dynamic pressure (q D ), 
which is related to the divergence velocity (U D ) by 



After evaluating the determinant, we obtain the following 
expression for q D . 


4D 


= ( 


66 ) 


r 

4^1^ 


16 


X 7 Z ll + I 7 I 8 ~ D„D 


11 66 


I 2 ] 
6 J 


2D 


[ I Vv - 


16 
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°H 3 6 


3a 


(2-43) 


When the integrals are evaluated using the beam mode 
shapes in Appendix B, the contribution of 

2 


D 11 C ~ 
48D 66* 


7 Vn 


will be small. The contribution of 




D U D 66 


to the value of q_, while not negligibly small, is much less 
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than that of 


2D 


16 


I 3 I 6 


'11 


The denominator of Eq. 2-44 will be zero when 


'16 

’ll 


eI 5 I 7 

2tl 3 I 6 


= 0 .3 — (for examples in 

this study) 


(2-44) 


When equals this critical value, q Q becomes infinite 

and the wing will not diverge. Above the critical value, 
q Q is negative and no real divergence occurs because U D will 
be an imaginary number. Below the critical value, including 
all negative values for D-^/D., 3 / there will be a finite, 
real U D . 

The physical mechanisms which cause this to be so are 
readily apparent in Fig. 3a and b. Figure 3a is a plate 
with a substantial positive D^. In this case, aerodynamic 
forces arising from plate displacements tend to return it 
to equilibrium, while just the opposite occurs with the plate 
having negative D lg (Fig. 3b). 

For the case where there is a known, nonzero, initial 
angle of attack, and a known dynamic pressure, the static 
deflections of the structure (a and w E ) can be obtained by 
solving the equations of Eq. 2-41 simultaneously. For this 
study, only tip deflections are evaluated, because they 
were obtained experimentally. 
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Eq. 2-41 can be rewritten as 


k ll q l + k 12 q 2 = a l a ( 


(2-45a) 


k 21 q l + k 22 q 2 " a 2 a ( 


(2-45b) 


where 


D ll c 


2 D... ac. 

16 „ l .. 

~V~ T 6 " qi 3a 1 3 
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, 3C. 
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qei 3a I 2 


By solving Eqs. 2-45a and b simultaneously, one obtains the 
following expressions for q^ and q^. 


a 2 a 0 ~ k 22 q 2 


(2-46a) 


[ a L k 21 ~ * 2 k ll 1 
[k 12 k 21 “ k ll k 22 J 


(2-46b) 


The generalized deflections are then transformed to 
actual deflections (a and w_) using Eqs. 2-9 and 2-10, 
where 4>^ and <$> a are evaluated at x = l. As before, a 
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(Fig. 2b) equals the sum of cXq and . The total angle of 
attack and bending deflection at the tip are 


<J> U=£) 

a TIP “ t 1 + ( 


a l k 21 ~ a 2 k ll 
k 12 k 21 “ k ll k 22 


) la. 


(2-47a) 


w T 


= [ 


'TIP 


2 _ _ 2_2 
k 21 *21 


a l k 21 ~ a 2 k ll 

1c V • V V 
12*21 *11 22 


) ] ^ (x=2.) a Q 
(2—4 7b) 


2.5 Free Vibration Problem 


The free vibration problem is formulated by setting 
the generalized forces (Q^) equal to zero in either Eqs. 
2-25a and b, or Eqs. 2-29a, b and c. To make the equations 
of motion into algebraic equations instead of differential 
equations, one assumes that the motions of the plate will 
be harmonic. For harmonic (sinusoidal) motion, the gener¬ 
alized displacements can be expressed as 


q = qe 


iwt 


« - 2 ,iwt 
q = -u qe 


(2-48) 


where co is the frequency. These expressions are substituted 
into the differential equations of motion (Eqs. 2-25 and 
2-29) to obtain the sinusoidal equations of motion. The 
sinusoidal equations of motion are presented here, in matrix 







form, for both the two term and three term deflection 
equations. 








40 


vibration frequencies (w n ) are obtained by taking the square 
root of each eigenvalue. From Eq. 2-49, one obtains two 
natural frequencies (first bending and first torsion), and 
from Eq. 2-50, one obtains three natural frequencies (first 
bending, second bending and first torsion). 

To simplify calculations, the warping stiffness term 



48 D 66* 2 


has been eliminated because its contribution is small for 
plates with moderate to high aspect ratios. The frequencies 
obtained using the mode shapes in Appendix B are compared 
in Chapter IV with experimentally obtained frequencies. 

2.6 Aeroelastic Flutter Problem 

The flutter problem is formulated using an indirect 
method widely known as the U-g method. In this method, 
the structural damping coefficient (g), introduced into 
the equations of motion, is plotted versus velocity for 
each vibration mode. Since solutions to the equations of 
motion represent conditions for neutral stability, the 
value of g obtained in this manner represents the amount 
of damping that must be added to the structure to attain 
neutral stability (flutter) at the given velocity. Therefore, 






negative values of structural damping indicate that the 
structure is stable. Flutter will occur when the artificial 
structural damping equals the actual damping of the structure 
To simplify calculations, the flutter problem will be formu¬ 
lated using only the equations of motion from the two term 
deflection equation. 

The equations of motion {Eqs. 2-24a and b) , as derived 
in Section 2.2, again assuming sinusoidal motion and neg¬ 
lecting the warping stiffness contribution, are 



\ V 3 *"* I 

[Q 2 /e iut | 
(2-51) 
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and the generalized forces (Q.) obtained from Eq. 2-20 are 
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Z 

Q 1 = / *i l e dx 

where L E 


c 

“2 


/ 2 p z d y 


(2-52a> 


l <p 

Q 2 " / -# «E dx 

where Mg 


_c 

c f 2 y p 2 d y 
2 


(2-52b) 


Since aerodynamic forces and moments obtained by assuming 
steady flow are generally inadequate for realistic flutter 
calculations, complete unsteady aerodynamic expressions 
will be used. Unsteady aerodynamic lift (Lg) and moment 
(Mg) expressions derived in Ref. 13 have been suitably trans 
formed to the coordinate system in Fig. 2. The reader is 
cautioned that the coordinate system used throughout this 
study differs from the standard aeroelastic coordinate 
system used in Ref. 13. The difference being the +w E here 
is -h in Ref. 13. The aerodyanmic expressions are 

Lg = 7Tpb 2 [-Wg + Ua - baa ] + 2;rpUbC (k) [-w £ 

+ Ua + b(1/2 - a) a ] (2-53a) 

Mg = Trpb 3 [~aWg - U(l/2 - a)a - b(l /8 + a 2 )a] + 

rrpUb 2 (1/2 + a) C(k) [-Wg + Ua + b(l/2 - a) a 1 


(2-53b) 
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where b 
a 

k 

C(k) 


semichord 

(distance EA is aft of midchord)/b (zero in 
this study) 

reduced frequency (wb/U) 

Theoaorsen function 


Assuming sinusoidal motion 


w E = w E e 
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a = e 


iwt 


w E = 


iwt 


a = lwae 
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w„ = -w 2 w_e 
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•• 2 a 

a = -w ae 
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the aerodynamic expressions, after some algebraic manipu¬ 
lation, become 
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[ M 3 + iM 4 ] = (1/8 + a 2 ) + - (1 - 2 ^ - C(k) 

+ i- [(1/2 - 2a 2 ) C (k) - (1/2 - a) ] 


Equations 2-54a and b are substituted into Eqs. 2-52a and b. 
One observes that, since the plate in this study has a con¬ 
stant chord, variables b, [ + il^ ), [ + iL^ ], 

[ + iM 2 ], and [ + iM 4 ] are not functions of x and 

can be brought outside the integral. Making use of the trans¬ 
formation equations (2-9 and 2-10), and expressing the integ¬ 
rals in nondimensional form, we arrive at the final expres¬ 
sions for the generalized forces. 
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Substituting Eqs. 2-56a and b into Eq. 2-51 and can¬ 
celing e lut yields a new form of the equations of motion, 
which are written here in contracted matrix form. 


[ K - od 2 A ] a * 0 


(2-56) 


where the elements of the K.. and A.. matrices are given as 
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Structural damping is proportional to, and opposes, 
the deflection of the structure. The structural damping 
coefficient (g) is introduced into Eq. 2-56 by multiplying 
the B matrix by (1 + ig). The methodology behind this step 
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is treated in depth in Ref. 13. An assumption here is that 
the linear damping coefficient (g^) is equal to the rotary 
damping coefficient (g^). It is also convenient to divide 
Eq. 2-56 by -w 2 and combine it with the (1 + ig) term to 
form a complex eigenvalue (Z), defined as 

Z 5 ( - + -- i - g ) (2-57) 

u > 2 

Equation 2-53 now becomes 

[ A - KZ ] q =0 (2-58) 

where the elements of g are real and the elements of 4 are 
complex. The matrix [ 4 - KZ 3 is put into standard eigen¬ 
value form by premultiplying by K -1 , resulting in 

[ K -1 A - IZ ] 2 = 0 (2-59) 

where I is an identity matrix. A unique solution to Eq. 

2-59 is obtained by setting the determinant of the coeffi¬ 
cient matrix equal to zero, and solving for Z. For a 2 x 2 
matrix, the two Z-values correspond to the two plate vibra¬ 
tion modes (first bending and first torsion). 

The solution procedure is to pick a value for the re¬ 
duced frequency (k) and solve Eq. 2-59 for Z^. The frequency 
structural damping coefficient and velocity are obtained 
from Z by the following relations. 
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u = —-— g = — { U * bw/k (2-60a,b,c) 

/EeTTF Re{ Z} 

A new k value is chosen and the procedure is repeated. 

Finally, a U-g diagram is created by plotting the structural 
damping for each root (Z) versus velocity, or a nondimensional 
flutter speed (U/bco^) , were ui a is the torsional natural vib¬ 
ration frequency. One can get a reasonable approximation for 

to by picking a very large value for k and solving Eqs. 2-59 
ct 

and 2-60a for the torsion root. This is essentially the 
structural vibration frequency when U equals zero. U-g 
diagrams for all the laminated plates considered in this 
study appear in Appendix B. The flutter velocity (Up) was 
conservatively chosen to be the point where the damping 
coefficient from either root first crosses the U-axis. 

Solving the flutter problem is further complicated by 
the presence of the Theodorsen function (C(k)). The Theo- 
dorsen function is a nonrational functional which can be 
expressed exactly in terms of modified Bessel functions of 
the second kind. Since Eq. 2-59 was to be implemented on 
a digital computer and solved many times for different 
values of k, an approximation for C(k) was desired. By 
taking the Laplace Transform of R.T. Jones' exponential 
approximation, one obtains a rational expression for C(k). 

This expression is 16 
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0.5s 2 + 0.2808s + 0.01365 

C {s) = -=- (2-61) 

+ 0.3455s + 0.01365 

where s is ik. 

The flutter problem is now completely defined. For 
actual computation, taking the determinant of the 2x2 
coefficient matrix is straightforward, the only complica¬ 
tions being that (1) the K-matrix must first be inverted 
and (2) that the coefficients of the A-matrix are all com¬ 
plex. Because of the simplicity of the problem, the temp¬ 
tation is to solve it manually. However, be advised that a 
typical U-g plot requires on the order of 20 data points. 
Each data point requires that Eq. 2-59 be solved for a new 
k value. 
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CHAPTER III 
EXPERIMENTS 


3.1 Test Specimen Selection 

Prior to making any test specimens, criteria which de¬ 
fined desirable, and sometimes essential, characteristics 
of the test specimens were established. These criteria in¬ 
cluded 

(1) Test specimens must exhibit a wide range of bending- 
torsion coupling stiffness. 

(2) Test specimens must be rectangular, constant thickness, 
flat plates with zero sweep. The requirement for a 
flat plate dictated using a symmetric laminate, since 
an unsymmetric laminate would warp during cure. 

(3) Test specimens would be made from Hercules AS/3501-6 
graphite/epoxy since it was available at M.I.T. 

(4) Test specimens would have an aspect ratio of four, 
since this would be representative of an aircraft wing 
with an aspect ratio of eight. 

(5) Test specimens should exhibit flutter and divergence 
within the 0-30 m/sec speed range of the acoustic 
wind tunnel at M.I.T. 

(6) The test specimens should be small enought to be made 
using standard TELAC (Technology Laboratory for Ad¬ 
vanced Composites) ply cutting templates and curing 
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plates. Additionally, a small size was desirable be¬ 
cause the effect of slight material warping, often 
encountered during curing, would be minimized. 

(7) The laminate should be able to withstand repeated 
large static and oscillatory loads. 

Criteria five and six directly oppose one another when 
trying to arrive at an optimum design. Trying to satisfy 
both of them dictated a very thin plate. Also, to achieve 
substantial bending-torsion coupling, an unbalanced lamina¬ 
tion scheme was chosen. Unbalanced laminates have an unequal 
number of +0 and -6 plies. However, a midplane symmetric 
ply arrangement was still required. Carrying the unbalance 
to the limit results in a unidrectional laminate with all 
plies at +0 or all plies at -9. To improve the toughness 
(criterion seven) of a unidirectional, off-axis laminate, the 
center two plies would have a zero-degree ply angle. The 
total thickness was chosen to be six plies. Five different 
laminates were selected as representative of a wide range 
of bending-torsion coupling stiffnesses, both positive and 
negative. The first two laminates, [+30-/0] and [-30-/0] , 

m S £ 9 

were chosen because a unidirectional laminate with 0 = 30° 
exhibits the highest bending-torsion coupling stiffness 
(D^g)^. The next two laminates, t-+45 ] s ana [-452/0] g , 

were chosen because a unidirectional laminate with 8 = 45° 
exhibits the highest torsion stiffness (D--), but still has 
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a large D^g. The final laminate, [ 0 ^/ 90 ] s , was chosen be¬ 
cause a unidirectional laminate with 0 = 0 ° exhibits the 
highest bending stiffness (D^) ar *d 2e ro D i 6 * During the 
test program, another laminate, [*45/0] , was added because 
it had the same theoretical and Dgg as [ + 452 / 0 ] , with 
less D^g. To minimize the number of test specimens which 
had to be constructed, the [+ 82 / 0 ] s and [-9 2 /0] s test spe¬ 
cimens were actually the same laminate, simply rotated 
180 degrees about the x-axis (Fig. 1) 

The test specimens would have to have an overall 
length of 330 mm (13 in) and a chord of 76 mm (3 in), be¬ 
cause this was the largest size that could be constructed 
satisfying criteria four and six. The overall length in¬ 
cluded a 25 mm (1 in) loading tab, making the effective 
cantilever plate length 305 mm (12 in). Finally, in an 
effort to minimize plate stiffness, an airfoil shaped 
fairing would not be used on the plate for wind tunnel 
tests. This had the added advantage of making the plate 
stiffness easier to calculate. 

3.2 Test Specimen Preparation 

The test specimens were constructed from Hercules 
AS/3501-6 graphite/epoxy prepreg tape from Lot No. 1643. 

The tape was 305 mm (12 in) wide and had a nominal thick¬ 
ness of 0.134 mm (0.00528 in). Individual plies were cut 
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to the proper size and angular orientation using aluminum 
templates, and assembled into 305 mm (12 in) by 356 mm 
(14 in) laminates. The laminates and curing materials (Ap¬ 
pendix C) were arranged on an aluminum curing plate especi¬ 
ally designed by students in TELAC for use with graphite/ 
epoxy laminates having a length and width as stated above. 

The laminates were cured in a Baron model BAC-35 autoclave 
using the cure cycle listed in Appendix C. After curing, 
the laminates were post-cured in a forced air circulation 
oven at 350° F for eight hours. After post-curing, a rec¬ 
tangular test specimen 330 mm (13 in) long and 76 mm (3 in) 
wide was cut from each of the four laminates using a diamond 
coated cutting wheel mounted on an automatic feed, milling 
machine. 

The thickness of a graphite/epoxy plate tended to 
vary over its surface. Therefore, the plate thickness was 
measured at several different locations and averaged. The 
same procedure was followed in measuring the length and 
width, although the variation was much less. The averaged 
measurements for each test specimen appear in Appendix D, 
along with the nominal values. The thickness variation 
between laminates (all were six ply) was 0.033 mm (0.0013 in) 
with the average being 0.807 mm (0.0318 in). This compared 
favorably to the nominal thickness of 0.804 mm (0.0317 in) 
for a six ply laminate. The laminates were also weighed 
on a triple beam balance, from which the material density 
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(p)for each test specimen was calculated. The average den- 

sity was 1.53 x 10 kg/m (0.0554 lb/in°) with a maximum 

variation of 0.06 x 10 3 kg/m 3 (0.0022 lb/in 3 ). This again 

3 3 

compared favorably with a nominal density of 1.52 x 10 kg/m 
(0.055 lb/in^) for graphite epoxy. 

Loading tabs 83 mm (3.25 in) x 25.5 mm (1 in) were 
machined from 2.4 mm (0.094 in) aluminum plate and bonded to 
the base of each test specimen with epoxy, cured at room 
temperature. The loading tabs were intended to aid in align¬ 
ing the test specimen in the clamping fixture and to prevent 
damage to the plate surface fibers. 

To get an indication of the lateral deflections the 
plate would experience during the wind tunnel tests, strain 
gauges were attached to the base of each test specimen at 
the midchord, as shown in Fig. 4. Two Micro Measurements 
EA-Q3-1873B-120 strain gauges, from Lot. No. R-A21AB02 with 
a gauge factor of 2.085, were attached to each test specimen 
(one on each side) to measure bending strain. Two BLH-SR4; 
FAED-25B-12-59; Serial No. 5-AE-SC strain gauge rosettes 
from Lot No. A-315 with a gauge factor of 2.01 were attached 
to each test specimen (one on each side) to measure torsion 
strain. The two bending gauges were wired together as a 
two-arm bridge circuit with three external lead wires approxi¬ 
mately 305 mm (12 in) long. The two torsion gauges were 
wired together as a four-arm bridge circuit with four external 
lead wires. Wiring the strain gauges on either side of the 
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plate together, in this manner, doubled the signal output 
for each channel (bending and torsion) and provided auto¬ 
matic temperature compensation. This temperature compen¬ 
sation was very important during the wind tunnel tests, where 
air flowing over an uncompensated gauge would cause a zero 
shift as the tunnel velocity was varied. The final step, 
after all solder connections had been made, was to coat the 
gauges with Micro Measurements M-Coat A, an air-drying poly¬ 
urethane, for protection. 

3.3 Test Apparatus and Procedure 
(a) Static Deflection Tests 

The static deflection test setup is shown in Fig. 5. 

It consisted of a 330 mm (13 in) x 508 mm (20 in) plywood 
base with six vertical steel rods approximately 760 mm 
(30 in) long. The test specimen was clamped in a vise 
machined from a 25.4 mm (1 in) x 152 mm (6 in) x 229 mm 
(9 in) aluminum block, which was bolted to the base of the 
test fixture. Two removable, low friction, pulleys were 
attached to the vertical rods such that a force or moment 
could be applied to the test specimen at any location along 
its length. Rulers, graduated in 32nd's of an inch, were 
also attached to the vertical rods to facilitate measuring 
the test specimen's edge deflections. A deflection indicator 
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was constructed out of balsa wood. It had needle pointers 
on the ends and cotton threads attached to its midpoint and 
ends. The threads, when routed over the pulleys and attached 
to weights, transferred a force or moment to the test speci¬ 
men (Fig. 5) . 

The deflection indicator was aligned with the lines 
scribed on the test specimen at the 3/4 length point, and 
the test specimen was clamped in the vise. The pulleys were 
clamped to the middle two vertical rods (Fig. 5a) at the 
proper height, and threads from the center of the deflec¬ 
tion indicator were routed over the pulleys. The rulers 
were adjusted to the proper height and zeroed with respect 
to the deflection indicator pointers. Weights of 10, 20, 

30, and 40 grams were successively attached to the threads, 
first to give positive deflections, then to give negative 
deflections. As each weight was attached, the readings 
from both pointers were recorded, along with the weight. 

Next, the pulleys were moved to the corner rods (Fig. 

5b) and the end threads were routed over them, so as to pro¬ 
vide a positive moment when the weights were attached. Weights 
of 10, 20, 33.9, and 43.9 grams were successively attached to 
each thread of the couple, and readings from both pointers 
were again recorded along with the weights. The pulleys were 
then switched to the corner rods on the opposite diagonal, 
and the procedure was repeated for a negative moment. 
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For each data point, the lateral deflection of the 
elastic axis (w £ ) and the rotation of the test specimen 
about the elastic axis (a) were calculated using the data 
reduction formulas in Appendix E. The lateral deflections 
obtained from the load test (Test 1) were plotted versus 
load, with the slope being the bending flexibility influence 
coefficient (c^) . The angular deflections (a) from Test 1 
were plotted versus load, yielding the bending-torsion 
coupling flexibility influence coefficient (c 2 ^)• The an¬ 
gular deflections obtained from the moment test (Test 2) 
were plotted versus moment, yielding the torsion flexibility 
influence coefficient (c 22 )• The lateral deflections ob¬ 
tained from Test 2 were plotted versus moment to obtain 
the other bending-tors ion coupling flexibility influence 
coefficient (c 12 ) • The flexibility influence coefficient 
plots for the four test specimens appear in Appendix E. 

The flexibility influence coefficients for a [+0 2 /O] s anc * 
a [-0 2 /O]g laminate will be the same, except for the sign 
on c 12 and c 2 ^, since they were physically the same test 
specimen. Finally, the flexibility influence coefficients 
for each test specimen were arranged in a 2 x 2 matrix (C) 
and inverted to obtain the stiffness influence coefficient 
matrix (K). The results of the static deflection tests are 
discussed in Section 3.4. 
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(b) Free Vibration Tests 

The free vibration test setup, shown in Fig. 6, used 
the same vise as the static deflection tests. The vice was 
suspended by four spring steel strips, which allowed it to 
translate along the z-axis (Fig. 2), but restricted motion 
in all other directions. The vise was rigidly attached to 
a horizontally mounted Ling model 420 shaker using a cylin¬ 
drical aluminum adapter. The shaker had a peak force of 
445 N (100 lbs) and a frequency range of approximately 5 to 
3000 Hz. Since the shaker was driven by an audio amplifier, 
there was a certain amount of distortion to the sine wave 
output signal below 20 Hz. For some vibration tests, an 
Endevco model 7701-50 "Isoshear" accelerometer was mounted 
to the vise by a threaded mounting stud, and an Endevco 
model 222B "Micro-miniature" accelerometer was mounted to 
the test specimen, within 25 mm of the base, using Eastman 
910 adhesive. The outputs of these accelerometers, after 
passing through special amplifiers, were connected to a 
Tektronix type 502 dual beam oscilloscope so that both out¬ 
puts could be displayed simultaneously. Finally, a digital 
signal counter was attached directly to the signal generator 
to provide an accurate frequency readout. 

The test specimen was aligned and clamped in the vise. 

A white paper screen was suspended behind the test specimen 
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Figure 6. 


Free Vibration Test Apparatus 
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to aid in observing the vibration modes. The test began 
with the signal generator set at, or below, 3 Hz. The fre¬ 
quency was slowly increased until the first bending mode was 
excited. The frequency was very finely adjusted until the 
vibration amplitude peaked, as observed visually, and the 
frequency was recorded. This procedure was repeated for 
the second bending and first torsion vibration modes. Since 
flutter is normally associated with the lower (bending and 
torsion) vibration modes, only the first bending, second 
bending and first torsion frequencies were obtained for each 
specimen. 

In an effort to improve the accuracy and repeatability 
of the frequencies, the dual beam oscilloscope was used to 
observe the accelerations of the vise and the test specimen. 
The resonance frequency was defined as the frequency where 
a 90° phase shift occurred between the two accelerometer 
traces. If a complete 90° phase shift did not occur, then 
the resonance frequency was taken to be the frequency where 
the maximum phase shift occurred. In an attempt to obtain 
a clear phase shift at the first torsion natural mode, the 
specimen-mounted accelerometer was moved to different chord- 
wise locations, all the way to the edge of the test specimen 
The results of the free vibration tests are discussed in 


Section 3.4. 
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(c) Wind Tunnel Tests 

All wind tunnel tests were performed in the M.I.T., 
Department of Aeronautics and Astronautics, acoustic wind 
tunnel. The acoustic wind tunnel is a continuous flow tun¬ 
nel with a 1.5 m (5 ft) x 2.3 m (7.5 ft) free jet test sec¬ 
tion 2.3 m (7.5 ft) long, located inside a large anechoic 
chamber. The tunnel was powered by a 100 HP motor giving 
it a continuously variable velocity range of 0 to 32 m/sec 
(0 to 105 ft/sec). The tunnel control panel was located 
outside the anechoic chamber. The velocity was controlled 
by two levers (coarse and fine speed control). The coarse 
lever controlled the motor field current and was variable 
in fixed step increments only. The fine lever controlled 
the motor shunt current and was continuously variable. 

The total current being drawn by the motor, along with 
voltage, were displayed at the control panel on analog 
gauges. There was a 400 amp current limitation on the motor, 
which was protected by circuit breakers. The tunnel velo¬ 
city was read from an alcohol manometer, calibrated in 
inches of alcohol, and located at the control panel. The 
manometer was connected to a pitot tube located slightly 
forward of the test section. A hydraulic motor, recessed 
into the floor of the test section, provided a means of 
mounting the test apparatus. 
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The test setup, shown in Fig. 7, consisted of a turn¬ 
table machined from aluminum mounted on a 635 mm (25 in) 
tall pedestal made of 152 mm (6 in) steel pipe. The pedes¬ 
tal was, in turn, mounted to the floor of the wind tunnel 
test section. The vise from the static deflection tests 
was fastened to the free rotating portion of the turntable. 

A wooden cover disk 508 mm (20 in) in diameter, shown re¬ 
moved in Fig. 7a and attached in Fig. 7b, was used to pro¬ 
vide smooth airflow past the test specimen. The disk, which 
had the angle of attack marked on its edge in 2° increments 
up to 18°, rotated with the test specimen. A pointer at¬ 
tached to the fixed base provided a consistent means of 
reading the angle of attack. An aluminum rod, attached to 
to the disk and extending outside the test section, allowed 
the angle of attack to be varied while the test was in 
progress. 

A terminal strip attached to the underside of the disk 
provided a convenient means of connecting the test specimen 
strain gauge lead wires to the instrumentation wiring harness. 
The instrumentation system was a Gould 2400 series four chan¬ 
nel strip chart recorder. The recorder had a four arm D.C. 
bridge preamplifier installed in channel 1 and a two arm 
D.C. bridge preamplifier installed in channel 2. The recorder 
also had an internal time base generator which marked the 
strip chart at one second intervals. The preamplifier had 
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the capability of both attenuating and amplifying the input 
signal. The strip chart speed could be set at 5, 10, 25, 50, 
100, and 200 mm/sec plus each speed could be divided by 100. 
The chart speed was accurate to * 0.5%. To provide a record 
of the test specimen deflections, the torsion strain gauge 
was attached to channel 1 and the bending strain gauge was 
attached to channel 2. 

The test specimen was aligned and clamped in the vise. 
The cover disk was installed, and all gaps were covered with 
tape. The strain gauge lead wires were attached to the 
terminal strip, and the strip chart recorder was adjusted 
according to the manufacturer's instructions. The test 
number and test specimen, along with test section tempera¬ 
tures and atmospheric pressure, were recorded on a data 
sheet. The first three natural vibration frequencies were 
recorded on the strip chart by giving the test specimen a 
bending deflection and letting it oscillate, then giving 
it a torsion deflection and letting it oscillate. Static 
calibrations were performed by giving the test specimen a 
series of known pure bending and torsion tip deflections 
and appropriately annotating the strip chart deflections. 

The angle of attack was adjusted to zero, and the wind 
tunnel was started and set at a low velocity (< 8 m/sec; 

25 ft/sec). After the velocity stabilized, a reading was 
taken from the alcohol manometer andannotated on the strip 
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chart, in inches of alchol. The angle of attach was then 
swept from 0° to 12° (sometimes to 18° if deemed appropri¬ 
ate) in 2° increments, while the velocity was held constant. 
Each angle of attack was annotated on the strip chart. The 
angle of attack was returned to zero and the tunnel speed 
increased slightly. The new manometer reading was annotated 
on the strip chart, and the angle of attack sweep was re¬ 
peated. The strip chart recorder was run at 1 mm/sec for 
most of the test. When oscillation was detected, a two 
second expanded record (100 or 200 mm/sec) was made. For 
some test specimens, at higher angles of attack, the oscil¬ 
lations would begin very abruptly. To get better data on 
the angle of attack, where flutter actually started, the 
angle of attack was increased in 1° increments in this 
region. This test procedure was repeated until either the 
the maximum tunnel speed was reached, or the test specimen 
deflections became excessive. At higher tunnel speeds, the 
maximum angle of attack tested was reduced to prevent damage 
to the test specimen. For the specimens which exhibited 
flutter at zero angle of attack, the normal procedure was 
to obtain two or three more data runs at velocities slightly 
higher than the velocity where the test specimen fluttered 
at zero angle of attack. 

Divergence testing often results in the destruction of 
a test specimen. Since this was not a desirable result. 
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wind tunnel speeds were increased in very small increments 
as the divergence velocity was approached. The deflections 
of the test specimen were closely monitored, and when it 
would no longer return to a zero deflection state, the 
tunnel velocity was not increased further. The extreme 
flexibility of the test specimens also increased the allow¬ 
able margin for error. 

Finally, after all the desired data was taken, the 
tunnel was shut down and the temperature of the test sec¬ 
tion was recorded on the data sheet. Since the wind tunnel 
was a continuous flow variety, temperature increases of up 
to 10° F were encountered during test runs. The test speci¬ 
men was removed from the vise and another test specimen 
installed. When going from a [+0 2 /°^ specimen to a 
[ — 0 2/0 3 s specimen, the turntable was simply rotated 180 
degrees, while the specimen remained undisturbed. The test 
procedure was then repeated for the new test specimen. 

After the wind tunnel tests were complete, the mano¬ 
meter readings were converted to indicated velocity, and 
then to true velocity using the conversions in Appendix E. 
The temperature used for the true velocity correction was 
the average of the start and finish temperature for each 
run. The oscillation frequencies were obtained manually 
and annotated on the strip chart, along with the indicated 
velocity. The static calibrations were used to obtain an 
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approximate relationship between test specimen deflection 
and strain gauge output using the procedure outlined in 
Appendix E. This relationship was used, on selected tests, 
to calculate the test static deflections (not oscillation 
deflections). The peak to peak oscillation amplitude of 
the bending and torsion strain gauges, on selected tests, 
was also annotated on the strip chart. Finally, the flutter 
boundary velocities, along with the dimensionless reduced 
flutter speeds (U/b<jJ a ) were plotted versus angle of attack 
for each test specimen. On selected tests, the test speci¬ 
men deflection, for a given angle of attack, was plotted 
versus velocity. Also, flutter amplitude, for a given 
angle of attack, was plotted versus velocity for selected 
tests. These test results are discussed in Section 3.4. 

3.4 last Results 

(a) Static Deflection Tests 

The flexibility influence coefficients (c^..) which are a 
measure of the structure's deflection for a given load, were 
obtained by determining the slope of the lines plotted in 
Appendix E. All the plots were linear, with the exception 
of c 2 2 f° r the [0 2 /90 ] s test specimen. For this plot, a 
linear approximation was used to determine the slope. The 
data scatter on all the plots was less than 10%. The flexi¬ 
bility influence coefficients could have been arranged in 
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a 2 x 2 matrix and inverted to obtain the stiffness influence 
coefficients (k^j)» which are a measure of the force a struc¬ 
ture exerts for a given displacement. However, for this 
section, the flexibility influence coefficients will be exa¬ 
mined directly. Table 1 lists the experimentally obtained 
inflence coefficients for each test specimen. 

TABLE 1 

EXPERIMENTAL FLEXIBILITY INFLUENCE COEFFICIENTS 

Test C 11 c 12 C 21 c 22 Coupling 

Specimen (ft/lb) (ft/ft-lb) (rad/lb) (rad/fflb) Factor (<) 


[0 2 /9° ] s 

0.18 

-0.03 

-0.02 

2.3 

0.0014 

[ ± 45/0] g 

0.60 

-0.16 

-0.13 

1.06 

0.033 

[ + 45 2 /°1 s 

0.98 

-0.76 

-0.75 

2.7 

0.22 

[- 45 2 /°] s 

0.98 

+ 0.76 

+ 0.75 

2.7 

0.22 

[+30 2 /0] s 

0.64 

-0.80 

-0.82 

3.0 

0.34 

[-3° 2 /°] s 

0.64 

+0.80 

+0.82 

3.0 

0.34 

A measure of 

the accuracy 

of these 

static deflection 


tests was that, according to Maxwell's Law of Reciprocal 
Deflections, c^ 2 must equal c 21 . All of the test specimens 
showed excellent agreement in this area. 

The coupling factor (tc) , included in Table 1, was de¬ 
fined as 

_ C 12 C 21 
< = - 

C 11 C 22 
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and was a dimensionless measure of the amount of bending- 
torsion coupling a test specimen possessed. One observes 
that going from a balanced to an unbalanced laminate sig¬ 
nificantly increased the coupling. Also, the [+30 2 /0] g and 
[-30 2 /0] s test specimens had the highest coupling of those 
tested. 

In going from the balanced [*45/0] test specimen to 

S 

the unbalanced [+45 2 /0] g test specimen, the torsional flexi¬ 
bility ( c 22 ) increased by a factor of 2.5, and the bending 
flexibility (c^) increased by a factor of 1.6. The reader 
should also note that the [+3O 2 /0] g test specimen was more 
flexible in torsion than the [O 2 /90] g test specimen. 

(b) Free Vibration Tests 

Excitation of the torsion vibration mode was extremely 
difficult on the [C> 2 /90] s anci [ ±4 5/0] g test specimens, 

which had very little bending-torsion coupling. Addition¬ 
ally, the torsion mode was only lightly excited on the other 
test specimens. The accelerometer bonded to the test speci¬ 
men was also largely unsuccessful for the torsion mode, be¬ 
cause the phase shift was very slight. Since the shaker 
was driven by an audio amplifier, the output signal was 
distorted below 20 Hz. This made using an oscilloscope to 
detect phase shift in the first bending mode difficult. 

The phase shift technique did, however, work well for the 








71 


second bending mode. 

Frequencies of the lowest three natural vibration modes, 
obtained from both the shaker and wind tunnel tests, are 
tabulated in Appendix E. The averaged frequencies for each 
test specimen appear in Table 2. Frequencies from both 
sources were all within 3% of the average. The 
and [—90] s test specimens, being the same physical plate, 
had the same natural frequencies. 


TABLE 2 

EXPERIMENTAL NATURAL FREQUENCIES (AVERAGE VALUES) 


Test 1st bending 2nd bending 1st torsion 

Specimen freq (Hz) freq (Hz) freq (Hz) 


[0 2 /9°] s 

11.1 

69.0 

42.2 

[*45/0] s 

6.07 

38.5 

77.4 

[+45 /0] ' 

) 



-VoJ 

> 4.83 

1 

30.1 

51.1 

[+30 /O], 




z s ( 

[-30-,/°]/ 

► 5.97 

1 

36.1 

58.4 


In going from the balanced [*45/01 test specimen to 

5 

the unbalanced [+45-/0] test specimen, the torsion fre- 

Cm S 

quency, which is largely dependent on the torsional flexi¬ 
bility of the specimen, decreased by 34%. Also, the first 
bending mode frequency decreased by 20%. These results 
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were consistent, at least in direction, with those of the 
static deflection tests. 

The situation becomes slightly more confusing when the 
frequencies of the other test specimens are compared with 
their flexibility influence coefficients. The first and 
second bending mode frequency trends correlate well with 
the trends of the bending flexibility influence coefficients 
However, the trends for the torsion frequencies and the tor¬ 
sion flexibility influence coefficients had poor correlation 
This indicated that the bending torsion coupling must be 
primarily affecting the torsion frequency. As an aside, 
when stiffness influence coefficients were compared to 
natural frequencies in this manner, one was led to conclude 
that the coupling primarily affected the bending frequencies 

All the test specimens, except the [O 2 /90] s specimen, 
exhibited first bending as the lowest frequency, followed by 
second bending, and finally first torsion. The [0 2 /90] s 
test specimen started with first bending, followed by first 
torsion and finally second bending. 

(c) Wind Tunnel Tests 

The primary objective of the wind tunnel tests was to 
determine the zero initial angle of attack (<*q = 0) diver¬ 
gence and flutter velocities for each test specimen. Since 
the test velocities were low and the test specimens quite 
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flexible, the onset of flutter at zero a Q tended to be gra¬ 
dual. For consistency, the flutter threshold was defined 
to be the velocity where the oscillation amplitude stabi¬ 
lized and the waveform became periodic. Determining the 
threshold for stall flutter was somewhat easier, because 
oscillation tended to begin abruptly when the stall angle 
of attack was reached. 

The type of flutter (bending, torsion or bending- 
torsion) could not be determined completely from the ampli¬ 
tude of the bending and torsion strip chart channels. This 
was because, for a plate with bending - torsion coupling, the 
elastic axis was not located at the midchord, where the 
strain gauges were located. So a pure bending deflection 
excited both the bending and torsion strain gauges. The 
problem was overcome by defining bending flutter as oscil¬ 
lations with a frequency very close to one of the bending 
natural frequencies. Torsion flutter was said to occur 
when the frequency was close to the torsion natural fre¬ 
quency, and bending-torsion flutter occured at a frequency 
between the other two. 

On test 9, with the 1-452/0] g test specimen, a torsion 
strain gauge failed due to excessive deflection about half 
way through the test. The test was completed using the 
bending channel only. The malfunction was in a wire sol¬ 
dered to the terminal, which was repaired and caused no 

further problems. A bending strain ga e on the (-30,/0] 

4 s 
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test specimen failed at the end of Test 3, also due to ex¬ 
cessive deflections. This gauge was not repairable, so 
Tests 14 and 15 were run with only one bending strain gauge. 
Temperature fluctuations during these two tests caused a 
zero drift on the bending channel. This zero drift made the 
data worthless for determining static deflections. However, 
the data were satisfactory for determining oscillation 
frequencies and amplitudes. 

The flutter threshold velocity was plotted versus ini¬ 
tial angle of attack in Fig. 8 . The test specimens exhibited 
markedly different behavior. The [0 2 /90] test specimen, with 
no bending-torsion coupling, exhibited bending-torsion flutter 
with a frequency of approximately 30 Hz for angles of attack 
below two degrees and torsion flutter with a frequency of 
approximately 42 Hz for higher angles of attack. The test 
specimens with large positive bending-torsion coupling, 

[-*-30 2 /° 1 3 anc * ( + 45 2 /0] s , exhibited primarily bending-torsion 
flutter at 28 Hz and 24 Hz, respectively. The flutter thres¬ 
hold velocity did not drop significantly with increasing 
angle of attack, probably because the coupling caused a de¬ 
crease in the tip angle of attack, preventing it from stall¬ 
ing. At an angle of attack of 18 degrees, the flutter thres¬ 
hold for the [+30-/0] test specimen dropped significantly 
and the flutter changed to torsion flutter at 55 Hz. The 
(±45 2 / 0 ] test specimen, which was stiff in torsion and had 
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only a small amount of positive bending-torsion coupling, 
would not flutter at zero angle of attack within the 32 m/sec 
maximum speed for the wind tunnel. However, at angles of 
attack above eight degrees, it behaved similarly to the 
[O2/90 ] s test specimen. The flutter was primarily torsion 
at 60 Hz. The test specimens with large negative bending- 
torsion coupling, [- 3 O 2 /O] and t“45 2 / 0 1 s , exhibited divergence 
at zero angle of attack. These two test specimens exhibited 
primarily bending flutter with frequencies of 4.5 Hz and 
4.0 Hz, respectively, at angles of attack greater than one 
degree. This flutter was more of a Von Karman vortex shed¬ 
ding buffet then a true stall flutter. 

Several data points in Fig. 8 have been marked with a 
(') to indicate a high frequency, low amplitude oscillation. 
The [O 2 /90 ] s test specimen exhibited bending oscillations 
at 69 Hz and torsion oscillations at 145 Hz. The [ ± 45_/0] 

4 S 

test specimen exhibited torsion oscillations at 108 Hz. 

These oscillations indicated that the second bending and 
second torsion modes had been excited, which happened because 
the plates were quite flexible. They do not represent the 
phenomenon with which this study is concerned. 

The effect of bending-torsion coupling on the zero 
angle of attack flutter velocity could not be determined 
directly from Fig. 8. This was because the test specimens 
had varying amounts of torsional stiffness, which strongly 













affect the flutter velocities. To cancel out the effect of 
torsional stiffness, the dimensionless reduced flutter 
velocity (U_/ba) ) was plotted versus angle of attack in 
Fig. 9. In this parameter, b is the semichord and is 
the first torsion natural frequency, which was obtained 
experimentally. 

The trends for each test specimen are similar to Fig. 

8 , as one would expect. However, in Fig. 9, the [0-/90] 

fc s 

test specimen had the highest zero angle of attack reduced 

flutter velocity, followed by the [+45 2 /0] s test specimen, 

and finally the [+30-/0] test specimen. Since the 

^ s 

[+30 2 /0] g test specimen had the largest positive bending- 
torsion coupling, the [+45 2 /0] g test specimen had somewhat 
less, and the [O 2 /90] g test specimen had none, one may 
conclude from this very limited sample size that zero angle 
of attack flutter reduced velocity was inversely propor¬ 
tional to positive bending-torsion coupling. 

The [ — 45 2 /0] s test specimen had a higher reduced di¬ 
vergence velocity (U D /bu a ) than the [ —30 2 /0] test specimen. 
So, using the same logic as before, one may conclude that 
divergence velocity was inversely proportional to increas¬ 
ingly negative bending-torsion coupling. Visual observa¬ 
tions indicated that divergence was impossible for test 
specimens with substantial positive bending-torsion 
coupling. 










Reduced Velocity ^ U/bw (dimensionless) 



Figure 9. Reduced Flutter Boundary Velocities 
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Static tip deflections, both a.. and w , for the 

fcl P e tip 

[+45-/0]. [*45/0] , and [-45-/0] were plotted versus velo- 
4 s s ^ s 

city in Figs. 10, 11 and 12, respectively. The tip deflec¬ 
tion data were obtained using the data reduction procedure 
in Appendix E. This data reduction procedure, which relates 
strain gauge readings to tip deflections through static cali¬ 
brations, was only an approximation. Therefore, Figs. 10 - 
12 must be examined for qualitative information only. The 
lines faired on the plots represent a "best fit" to the 
available data, which have been arbitrarily extrapolated to 
indicate some trends. The ot. . was the sum of the initial 

tip 

angle of attack (a-) and the elastic twist (a ). 

0 e 

Figure 10 indicates a linear relationship between 
lateral tip deflection and velocity. This seems odd, be¬ 
cause the force on the wing is a function of the dynamic 
pressure (q), which is a function of the square of the 
velocity. One would expect this relationship to be para¬ 
bolic. Two possible explanations are offered. First, the 
relationship was parabolic and the data were somewhat in 
error. Second, the bending-torsion coupling reduced the 
tip angle of attack to a low value as the test specimen was 
bent laterally. Unfortunately, the plot of the tip angles 
of attack doesn't support this second explanation. However, 
values for obtained from the strip chart were extremely 










Figure 10. Experimental Static Aeroelastic Deflection, [+45VO] Te 
Specimen s 
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Figure 11. Experirrental Static Aeroelastic Deflections, [ ± 45/0] s Itest 
Specimen 
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sensitive to slight errors in reading the strip chart. This 
was compounded by the fact that buffeting of the test speci¬ 
mens made it difficult to get an accurate strip chart reading. 

Figure 11 indicates a more parabolic relationship be¬ 
tween the lateral tip deflections and velocity. Also the tip 
angle of attack increases with increasing velocity as one 
would expect for a wing with low bending-torsion coupling. 

More data points were available for Fig. 11 and they have 
good continuity. 

Figure 12 shows the expected divergent response ten¬ 
dencies of a wing with negative bending-torsion coupling. 

The flutter amplitude was plotted versus velocity for 
the [+45-/01 , [*45/0] , and [-45-/0] test specimens in 

S S c, S 

Figs. 13, 14, and 15, respectively. The flutter amplitude 
was unreduced, unamplified, peak to peak oscillation data 
directly from the strip chart. It has no significance 
other than to indicate where the flutter starts, for a given 
initial angle of attack, and to indicate how rapidly the 
magnitude increases. For both the [ + 45 2 /0J s and £ — 45 2 /0] s 
test specimens flutter beings quite close together for 
several angles of attack. While, for the [ ± 45/0] g test 
specimen the occurrence is more spread out. 
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CHAPTER IV 

COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 

18 19 

Work done by Boyce , Crawley , and, in particular, 
Turner 10 showed that stiffness properties for graphite/epoxy 
laminates experimentally obtained from in-plane extensional 
loading were generally higher than material properties ob¬ 
tained from out-of-plane flexural loading. The standard 
TELAC values for E^, E^, G LT and v^, which appear in Ap¬ 
pendix A, were obtained for Hercules AS/3501-6 G/E from in¬ 
plane loads. Turner, through a series of out-of-plane 
flexural static and dynamic tests, established an alterna¬ 
tive set of material property values for AS/3501-6 G/E. 

These material property values also appear in Appendix A. 
Turner’s values, while not statistically significant like 
the TELAC values, do constitute a basis for comparison. 
Therefore, in this chapter, experimental results will be 
compared to theoretical results using both TELAC and Turner 
material properties for graphite/epoxy. Additionally, dis¬ 
cussion, in several instances, will be limited to quali¬ 
tative "trends" due primarily to the approximate nature of 
the data being compared. 

4.1 Static Deflection Tests 


The theoretical and experimental flexibility influence 
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coefficients for each test specimen appear in Table 3. For 
clarity, they have been arranged in 2 x 2 matrices. As 
stated before, c^ is the bending flexibility, c^ 2 is the 
torsion flexibility, and c^ 2 , are the bending torsion 

coupling flexibilities. 

Values for c^ 1 calculated from Turner material proper¬ 
ties were within 10% of experimental values. The c^'s 
from TELAC material properties were lower than Turner 
and experimental c-^'s, indicating that TELAC material pro¬ 
perties resulted in a plate that was too stiff in pure 
bending. The coupling flexibility influence coefficients 
obtained from TELAC data were, however, in better agreement 
with experimental data than those from Turner data. 

The torsion influence coefficient for the [*45/0] test 

s 

specimen obtained from Turner data had reasonable agreement 
with the experimental value. The torsion influence coef¬ 
ficients obtained from both Turner and TELAC data had poor 
agreement with experimental values for test specimens with 
large bending-torsion coupling. With one notable exception, 
all theoretical values for c 22 were lower than experimental 
values. The exception was the [O 2 /90] test specimen, where 
theory predicted a substantially higher value for c 22 than 
was obtained from the experiments. For all test specimens, 
except [O 2 /90] s , torsion influence coefficients obtained 
using Turner data had better agreement with experimental 
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values than did torsion influence coefficients obtained 
using TELAC data. 

There were several possible reasons for differences 
between theoretical and experimental torsion influence co¬ 
efficients. First, the two term deflection equation may 
not adequately represent the actual deflection of the 
structure. Additional terms would allow more degrees of 
freedom. Second, the assumption that the plate was chord- 
wise rigid when undergoing torsional deflections may not 
be true, especially for test specimens with large bending- 
torsion coupling. Third, the theory used to obtain the 
flexural moduli (D.^) for the plates, as presented in 
Chapter II, assumed small deflections and, consequently, 
neglected transverse shear. If transverse shear was pre¬ 
sent, the plate would become more flexible in torsion. 
Further, if transverse shear increased with increasing 
bending-torsion coupling, then for test specimens with 
large bending-torsion coupling, one would expect the experi¬ 
mental torsion flexibility influence coefficients to be 
much larger than the theoretical influence coefficients. 

This was the case in this study. 

The final significant trend was that the coupling 
factor (k) for the theoretical flexibility influence co¬ 
efficients was substantially higher than the experimental 
coupling factor for all test specimens except [Oj/SOlg. 
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Both the theoretical and experimental coupling factors were 
essentially zero for the [C>2/90] s test specimen. 


4.2 Free Vibration Tests 


A summary of the theoretical and experimental natural 

frequencies for the first three natural vibration inodes of 

each test specimen appears in Table 4. Also included in 

Table 4 are the theoretical uncoupled natural frequencies 

19 

calculated using the following equations. 


(first bending) 


(second bending) 


(first torsion) 


“IB 



“2B 



“IT 


1.57 


48D 66 

m£ 2 c 2 


(4-la) 

(4-lb) 

(4-lc) 


where is the bending flexural modulus, Dgg is the tor¬ 
sion flexural modulus, m is the mass per area, l is the 
length, c is the chord, and oj is the frequency in rad/sec. 

The theoretical first bending natural frequencies, 
calculated using Turner material properties, were all within 
5% of the experimental values, which was good agreement. 

The first bending frequencies calculated from TELAC material 
properties had poorer agreement with experimental values. 
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In all cases, the TELAC first bending frequencies were 
higher than experimental frequencies. The theoretical un¬ 
coupled first bending frequencies, as one would expect, 
showed good correlation with experimental frequencies for 
test specimens with little bending-torsion coupling, and 
poor correlation for test specimens with large bending- 
torsion coupling. 

The theoretical values for the second bending fre¬ 
quencies from Turner material properties showed reasonable 
agreement with the experimental values for test specimens 
with low coupling, and poor agreement for test specimens 
with large coupling. Second bending frequencies calculated 
from TELAC material properties were not as good as the 
frequencies from Turner data. The theoretical uncoupled 
values for the second bending frequency had generally 
poorer agreement with the experimental values than the 
coupled solutions did. 

The theoretical first torsion natural frequencies had 
generally poor agreement with experimental values. Errors 
for first torsion natural frequencies calculated using 
Turner material properties were lowest for the [*45/0] g test 
specimen and highest for the [+45 2 /0] s test specimen. More 
importantly, one observes that the theoretical torsion fre¬ 
quency values for the [^45/0] and [+45-/0] test specimens 
were essentially the same, while the experimental torsion 













frequency for the t+45 2 /0] g test specimens was substantially 
lower than the experimental torsion frequency for the [*45/0] 
test specimen. The higher torsion frequency for the [*45/0] 
test specimen indicated that it was much stiffer in torsion 
than the [+45-/0] test specimen. As before, this was pos- 

Z S 

sibly due to neglecting transverse shear in the problem 
derivation. Also, the fact that only the first three vib¬ 
ration mode shapes were used to formulate the assumed deflec¬ 
tion equation could affect the solutions somewhat. Finally, 
as before, the assumption that the plate was chordwise rigid 
during torsion vibrations could also affect the solution. 

The theoretical natural frequencies should possibly be 
checked by a more extensive Rayleigh-Ritz formulation and 
a finite element analysis. 

4.3 Divergence Velocities 

Theoretical and experimental divergence velocities for 
each test specimen appear in Table 5. Velocities were ob¬ 
tained for plates using both TELAC and Turner material prop¬ 
erties, and for both a two dimensional lift curve slope 
(2c^/2a) and one empirically corrected for finite span using 
the correction expression in Chapter II. Of the six speci¬ 
mens tested, only three would possibly exhibit any divergence 
at zero initial angle of attack. These three test specimens 
were [O 2 /90] s , [ — 45 2 /0] g , and [-30 2 /0] g . 


Of these three. 
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the [0-/90] test specimen exhibited flutter prior to diver- 

^ S 

gence, so no experimental divergence characteristics were 
obtained. However, a theoretical divergence velocity was 
still calculated for that specimen. 

Theoretical divergence velocities were lower than experi¬ 
mental velocities in all cases. The finite span divergence 
for the [-45 2 /0 ] s and [-30 2 /0] g test specimens, calculated 
using Turner material properties, were lower than the diver¬ 
gence velocities calculated using TELAC material properties. 

An experimental divergence velocity was not obtained for 
the [O 2 /90] g test specimen. However, the wind tunnel was 
taken up to 27 m/sec (90 ft/sec) with this test specimen 
and divergence did not occur. This contrasts to a theoreti¬ 
cal finite span divergence velocity, calculated with Turner 
material properties, of 22 m/sec (71 ft/sec). This indi¬ 
cated that the theory would not adequately model the [O 2 /90] g 
test specimen. 

The fact that the theory predicted a lower divergence 
velocity than found from experiments would suggest that the 
actual plates were stiffer than the theoretical plates. This 
observation is in sharp contrast to the observations in Sec¬ 
tions 4.1 and 4.2, where the theoretical plates appeared 
stiffer. The observation would tend to rebut the explana¬ 
tions given in those sections. Since divergence is a func¬ 
tion of aerodynamic force on the plate, as well as plate 
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stiffness, the aerodynamic strip theory needed to calculate 
the loads on the plate may have been a source of error. 
Further, the finite span correction may have resulted in the 
theoretical air loads being greater than the actual air loads. 
Reference 13 suggests that an alternate empirical correction 
for finite span effects may provide better correlation be¬ 
tween theoretical and experimental results for divergence 
problemr. That empirical correction is 

—^ a a- E ——— ] (4-2) 

3a AR + 4 

Using this finite span correction, theoretical divergence 
velocities were raised slightly. A final source of error 
may have been the criteria used during the wind tunnel 
tests to determine the divergence velocity. This, however, 
is an improbable source of error. 

4.4 Static Aeroelastic Tip Deflections 

The theoretical tip lateral deflections (w-, ) and 

t-Tip 

tip angles of attack (a TI p) were calculated using Turner 
material properties, the empirical finite span correction, 
and the solution technique in Chapter II. The tip deflec¬ 
tions were plotted versus velocity for the [+45-/0] 

S 

[-45/0] g , and [-45 2 /0] g test specimens in Figs. 16, 17, and 
18, respectively. To allow a comparison with experimental 
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results, the approximate experimental tip deflections, dis¬ 
cussed in Chapter III, were included in these figures. The 
lines represent the theory and the plotting symbols are the 
experimental values. Since data from both theory and experi¬ 
ment were approximate, a valid comparison must be limited 
to trends only. 

In Fig. 16, the lateral deflection trends for the 

[+45-/0] test specimen show reasonable agreement between 
z s 

theory and experiment, especially at higher velocities. 

The trends for the theoretical and experimental tip angle 
of attack, however, did not agree. The experimental values 
indicated a very slight increase in tip angle of attack as 
velocity was increased. Theory, on the other hand, showed 
a substantial decrease in tip angle of attack as velocity 
was increased. The implications of the theoretical results 
were that the torsional moment resulting from bending-torsion 
coupling more than overcame the aerodynamic moment. A pos¬ 
sible reason was that the theoretical coupling factor (k), 
discussed in Section 4,1, was much larger than the experi¬ 
mental k. Of course, there was also the problem with read¬ 
ing the strip chart, as discussed in Chapter III. 

The agreement between theoretical and experimental 
lateral tip deflections for the [*45/0] test specimen in 

S 

Fig, 17 was not quite as good as for the [ +-45 2 /° 1 s test 
specimen, but was still reasonable. Again, the trends for 


1 












2 


the tip angles of attack did not correlate between theory 

and experiment. The theoretical tip angles of attack did 

not decrease as much for the [*45/0] test specimen as they 

did for the [+45-/0] test specimen, but the experimental 

fc s 

tip angle of attack values increased more for the [*45/0] 

s 

test specimen. In view of the fact that the [*45/0] test 
specimen possessed slight positive bending-torsion coupling, 
one would anticipate this type of behavior. The reasons 
for the lack of correlation between experiment and theory 
then would be the same as for the [+452/0] g test specimen. 

Divergent tendencies of both the tip lateral deflec¬ 
tions and angles of attack for the [-45 2 /0] g te sh specimen 
are apparent in Fig. 18. The trends were similar for both 
the theoretical and experimental results, although the ex¬ 
perimental data did not correlate well with theory for indi¬ 
vidual angles of attack. Finally, the experimental tip angle 
of attack data appeared to diverge faster than the theore¬ 
tical data. 

4.5 Flutter Velocities 

The theoretical and experimental flutter velocities at 
zero angle of attack for each test specimen appear in Table 
6. Additionally, the nondimensional reduced flutter velo¬ 
cities for each test specimen are presented. Since the 
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TABLE 6 

a. FLUTTER VELOCITIES 


Test 

Specimen 

Theoretical 

TELAC 

(U-g Analysis) 
Turne r 

Experimental 

" 

[o 2 / 9 °] s 

18,0 

(59,0 

m/sec 
ft/secl 

17,4 m/sec 
(57.0 ft/sec) 

25.9 m/sec 
(85.0 ft/sec) 

[±45/0] s 

42.1 

(138) 


36.6 

(120) 

32.0 

(>105) 

[ + 45 2 /°] s 

40.2 

(132) 


35.1 

(115) 

27.7 

(91.0) 

t-45 2 /0] s 

50.3 

(165) 


43.6 

(143) 

diverged 

first 

[ + 30 2 /0] s 

34.7 

(114) 


30.5 

(100) 

27.1 

(89.0) 

[-30 3 /0] s 

50,0 

(164) 


42.3 

(142) 

diverged 

first 


b. REDUCED FLUTTER VELOCITIES 


Test 

Specimen 

TELAC 

Theoretical 

Turner 

Experimental 

[o 2 / 9 °] s 

2.21 


2.21 

2.56 

[*45/0] g 

2.21 


2.21 

>1.6 

[ + 45 2 /°] s 

2.11 


2.11 

2.24 

[-45^/0] 
l s 

2.64 


2.63 

diverged 
first (1.01) 

[+ 30 2 /01 s 

2.04 


2.04 

1.94 


C-30 2 /0] a 


2.93 


2.90 


diverged 
first (.839) 
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[- 452 /O] and I-302/0J s test specimens diverged prior to 
fluttering, no experimental flutter velocities were obtained 
for them. However, theoretical flutter velocities were 
calculated for these two test specimens. Also, the [ ± 45/0] s 
test specimen did not flutter within the speed range of the 
wind tunnel. 

The correlation between theoretical and experimental 
flutter velocities was generally poor. The flutter velocity 
for the [O2/90] g test specimen, calculated using Turner 
material properties, was much lower than the experimental 
flutter velocity. On the other hand, the theoretical flutter 
velocities for the [+45 2 /0] g and [+30 2 /0] g test specimens, 
calculated using Turner material properties, were much 
higher than the experimental flutter velocities. The flutter 
velocities calculated using TELAC material properties had 
generally worse correlation. 

The errors in the flutter calculations appeared to be 
directly related to the inability of the theory to accurately 
determine the torsion vibration frequencies. Therefore, one 
would expect a somewhat better correlation between theoreti¬ 
cal and experimental reduced flutter velocities (Up/bw^), 
since the contribution of the torsion vibration frequency 
had been eliminated. In fact, that was exactly what hap¬ 
pened. The [0 2 /90] s test specimen's error was reduced sig¬ 
nificantly, and the errors associated with the [+45-/0] 

2 s 
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and [+30-/0] test specimens were less than 10%. The 
z s 

reader should note that the ui used for the theoretical 
calculations was obtained from the U-g analsyis for a very 
large reduced frequency (k) value. 

Since the trends were reasonably well predicted by the 
theory, some conclusions on how bending-torsion coupling 
affects flutter velocity were drawn, based on theoretical 
results. The first conclusion was that increasingly posi¬ 
tive bending-torsion coupling reduced the flutter velocity. 
The second conclusion was that increasingly negative bending- 
torsion coupling raised the flutter velocity. This was an 
unrealistic situation, however, because the test specimens 
with negative coupling diverged long before they reached 
the flutter velocity. 
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CHAPTER V 

CONCLUSIONS AND RECOMMENDATIONS 

5.1 Conclusions 

The present study has sought to investigate the flutter 
and divergence behavior of graphite/epoxy plates with various 
amounts of bending-torsion stiffness coupling. As a conse¬ 
quence, the static deflection and vibration behavior has also 
been examined. 

There was reasonable agreement between theory and exp¬ 
eriments for the bending stiffness and first bending mode 
natural frequencies. Also, for test specimens with low 
bending-torsion coupling, there was reasonable agreement 
for the torsional stiffness and first torsion mode natural 
frequencies. There was, however, poor agreement between 
theory and experiments for torsion stiffness and frequencies 
of test specimens with large bending-torsion coupling. The 
poor agreement was probably due to an inadequate number of 
modes in the assumed deflection equation. Other possible 
sources of error were the assumption that the plate was 
chordwise rigid and the neglecting of the transverse shear 
effect. 

In this study, only the plates with negative bending- 
torsion coupling (-D^g) exhibited divergence within the 
speed range of the wind tunnel. For these test specimens, 
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the theory had reasonable agreement with the experiments. 

The errors in the theory were possibly due to the approxi¬ 
mate aerodynamic strip theory. For test specimens with 
positive bending-torsion coupling, the theory predicted that 
divergence would not occur. Experimentally, they did not 
exhibit a divergence tendency within the speed range of the 
wind tunnel. Based on the theory and the limited set of 
experiments, the following observations can be made. The 
zero angle of attack divergence velocity, for the plates in 
this study, became infinite for positive bending-torsion 
coupling above a certain critical value. This divergence 
velocity decreased as the bending-torsion coupling decreased 
from the critical value to zero, and continued to decrease 
as the bending-torsion coupling became increasingly negative. 

The theory and experiments had poor agreement for zero 
angle of attack flutter velocities. However, the agreement 
between theory and experiments for reduced flutter velo¬ 
cities (U/bw^) was good. This indicated that the inaccu¬ 
racies in the theory were in the prediction of the torsion 
vibration frequency. The test specimens with negative 
bending-torsion coupling diverged prior to reaching their 
flutter velocities, but a theoretical flutter velocity was 
still calculated for them. As before, based on the theory 
and the limited set of experiments, the following observa¬ 
tions can be made. 
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The zero angle of attack flutter velocity for the plate in 
this study was low for positive bending-torsion coupling. 

The flutter velocity increased as the positive bending- 
torsion coupling decreased to zero, and continued to in¬ 
crease as the bending-torsion coupling became increasingly 
negative. This latter situation, for negative bending- 
torsion coupling, becomes academic if divergence occurs 
prior to flutter. 

The test specimens exhibited markedly different flutter 
behavior as angle of attack was increased. The test speci¬ 
mens with low bending-torsion coupling exhibited classical 
bending-torsion flutter at zero angle of attack, but it 
changed to pure torsion flutter as the angle of attack was 
increased. Also, the flutter boundary velocity decreased 
rapidly with increasing angle of attack, stabilizing around 

g 

a = 12°, This trend is similar to the results Rainey ob¬ 
tained for metal plates and is classical stall flutter. The 
test specimens with large positive bending-torsion coupling 
exhibited bending-torsion flutter throughout most of the 
angle of attack range. Additionally, the flutter boundary 
velocities did not decrease with increasing angle of attack 
below a = 12°. This behavior was probably due to the fact 
that the tip of the plate tended to have a lower angle of 
attack than the root, so it did not stall. The test specimens 











109 

with large negative bending-torsion coupling diverged at 
zero angle of attack and exhibited a low frequency bending 
flutter for angles of attack greater than zero. The flutter 
boundary velocities here were significantly below the velo¬ 
cities for the other test specimens for angles of attack 
less than 12°. These bending flutter velocities decreased 
moderately with increasing angle of attack. 

5.2 Recommendations 

The theory should be further developed to give better 
results for the torsion vibration frequencies and flutter 
velocities. Improvement could probably be most easily faci¬ 
litated by going to a finite element analysis which allows 
more degrees of freedom and accounts for transverse shear. 

If a Rayleigh-Ritz formulation is used, more terms should 
be included in the deflection equation. For divergence 
velocity calculations, a better aerodynamic theory, such 
as lifting line theory, should be used. Also, the stall 
characteristics of the plate should be explored further. 

The effect of bending-torsion coupling as a passive 
gust alleviation technique might hold some promise. Further, 
bending-torsion coupling might be useful for increasing the 
aileron reversal speed for a swept wing. Finally, the effect 
of bending-torsion coupling on an oblique wing, particularly 
with the roll degree of freedom, might prove interesting. 
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APPENDIX A 

MATERIAL PROPERTIES 


TABLE A-l 

ORTHOTROPIC ENGINEERING CONSTANTS 
Hercules AS/3501-6 Graphite/Epoxy 


Property 


TELAC Values Turner Values 

Cin-plane loads) (out-of-plane loads) 





°LT 


Ply thickness 
Density 


130 x 10 9 Pa 
10.5 x 10 9 
0.28 

6.00 x 10 9 

1.34 x 10 -6 m 
1.52 kg/m 9 


98 x 10 9 Pa 
7.9 x 10 9 
0.28 

5.6 x 10 9 
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TABLE A-2 

LAMINATE FLEXURAL MODULI 


Laminate Modulus TELAC Turner 

Material Properties Material Properties 



D 11 

5.474 Nm 

“ ‘ . 1 

4.038 lb•ft 

4.126 Nm 

3.043 lb•ft 

[0 2 / 9 °] s 

°16 

0 

0 

0 

0 


°66 

0.2600 

0.1918 

0.2425 

0.1789 


°11 

1,996 

1.472 

1.550 

1.143 

[*45/0] 

s 

°16 

0.5789 

0.427 

0.4364 

0.3219 


°66 

1.422 

1.049 

1.074 

0.7920 


°11 

1.996 

1.472 

1.550 

1.143 

t + 4 V 0] s 

°16 

1.254 

0.925 

0.9456 

0.6975 


°66 

1.422 

1.049 

1.074 

0.7920 


D 11 

3.541 

2.612 

2.703 

1.994 

[+30 2 /0] s 

°16 

1,589 

1.172 

1.179 

0.8695 


D 66 

1.132 

0.8350 

0.8662 

0.6389 
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B.l Mode 
* x (x) = 

where £ ^ = 
a l 

<p 2 (x) = 

where e 2 
a 2 

<t> 1 (Z) = 

<P 2 (z) = 

»a ( *> - 


APPENDIX B 
THEORETICAL ANALYSIS 


Shapes 


14,15 


cosh £l (j) - cos £;l (j) - ai [sinh £l (2) - sine 1 (j)] 

(B-l) 


1.8751 

0,7341 


coshe 2 (^) - cose 2 (j) - a 2 I sinh e 2 { f ) “ sine 2 (f)l 

(B-2) 

4.6941 

1.0185 

<J> Cx) = sin(^y) (B-3) 

21 


2.00 

4> 1 (0.75fi.) 

= 1.32 

2.00 

<j> 2 (0.75£,) 

= 0.269 

.00 

<P (0.75 Z) 
r a 

= 0.924 


The following integrals were evaluated numerically using 
ten point trapezoidal integration. 








1 r l * a 

7 J 4>-,dx 

* 0 1 


1 r* * j 

1 / *a dx 

0 


= 0.783 


= 0.637 


J 3 = i / * 1 *« dx = °- 678 


I 4 H I 0 /A *1 dx = 1 * 00 


I ^ *a dx = 0* 500 


£2 o f *£*a dx * 3 ’ 76 


£3 g 7 (< t >”) 2 dx = 12.4 


1 0 f (<^) 2 dx = 1.23 


12 Q f *2*o dx = _6 ’ 43 


£3 O^ (< ^2 5 2 dx = 485 


11 = 


'12 = 


</ (<*>") dx = 3,04 


t </ *2 d * * i - 00 















:ural Danping ^ g (dimensionless) Structural Danping ^ g (dimensionless) 


w 


« . 


Figure B-l. U-g Diagram, [O^SOJg Test Specimen 





















Figure B-2. U-g Diagram, [*45/0] Test Specimen 











Structural Danping ^ g (dinensionless) Structural Danping ^ g (dimensionless) 
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Figure B-4. U-g biagram, [+30^/0] g Test Specimen 

















f — r 












Structural Danping ^ g (dimensionless) Structural Danping ^ o (dimensionless) 

















APPENDIX C 

GRAPHITE/EPOXY CURING TECHNIQUES 


3 Vacuum bag 

Air bleeder 
a (fiberglass cloth) 

3 Porous Teflon 

1 Aluminum plate 

J (top) 
a Nonporous Teflon 

Paper bleeder 

3 Porous Teflon 


6 -ply G/E laminate 
(with peel ply) 

Porous Tteflon 

Paper bleeder 

Nonporous Teflon 

Aluminum plate 
(bottom) 


Figure C-l. Curing Materials 
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TABLE C-l 

AUTOCLAVE CURE CYCLE 


1. Place aluminum curing plate with curing materials (Fig. 
C-l) in autoclave and attach vacuum hose. 

2. Perform vacuum tests Apply full vacuum (30 in) to curing 
plate, turn off vacuum pump, and watch bleed-down for 
five minutes. 

3. Seal autoclave, apply full vacuum, and raise autoclave 
pressure to 85 psi. 

4. Raise autoclave temperature to 240° F. 

5. Maintain 85 psi and 240° F for one hour. 

6 . Raise pressure to 90 psi. 

7. Raise temperature to 35Q° F. 

8 . Maintain 90 psi and 350° F for two hours. 

9. Decrease temperature at a rate not exceeding 5° per 

minute until temperature reaches 130° F. 

10. Release pressure and vacuum. 

11. Open autoclave and remove curing plate. 
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APPENDIX D 

TEST SPECIMEN DIMENSIONS 

TABLE D-l 

NOMINAL TEST SPECIMEN DIMENSIONS 


Parameter Value 

Length 12,0 in 304.8 mm 

Width 3.00 76.2 

Ply thickness 0.00528 0.134 

(nominal TELAC value) 

Total thickness 0.0317 0.804 

(6 x ply thickness) 


Density 0.055 lb/in 3 1.52 x 10 3 kg/m 3 

(nominal TELAC value) 










EXPERIMENTAL TEST SPECIMEN AVERAGE MEASUREMENTS 
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APPENDIX E 

TEST DATA AND DATA REDUCTION METHODS 

E.1 Flexibility Influence Coefficient Data Reduction 

Test apparatus characteristics: 

moment arm (1/2 couple) = 3.5 in 

distance between front and rear rulers = 8.0 in 

Test 1 raw data: 

applied forces (F^): ± 0,022 lb, ± 0.044, ± 0.066, 

± Q.088 

obtain front and rear deflections for each applied force 
Test 2 raw data: 

applied moments (M^): ± 0,013 ft*lb, ± 0.026, ± 0.044, 

± 0.056 

obtain front and rear deflections for each applied moment 
Data reduction for each test point: 

front + rear 

lateral deflection (w„) in units of feet = [ - 

E 768 

angular deflection (a) in units of radians = 

front - rear 


tan ^ [ 


256 


] 
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w T 


'11 


'21 


a W E 

; C 12 * ; C 22 


a 

Ml 


Plots of all test points for each test specimen appear in 
Figs. E-l to E-4. 

E.2 Indicated and True Velocity Calculations 
Raw data: 

manometer reading (h al(= ) in inches of alcohol 
temperature (°F) at start and finish 
barometric pressure (h^) in inches of mercury 

Indicated velocity calculation: 

specific gravity (Sg) of alcohol = 0.805 at 80° F 
specific weight (u) alc ) Of alcohol = Sg(a) H Q ) = 

(0,805) (62.4 lb/ft 3 ) = 50.2 lb/ft 3 


U. .. 
indicated 


2ui alc h alc 

P SL 


= 59.3 


h alc (inches) 


True velocity calculation: 


“true “ “indicated 1 > " here ® * 

P SL 1 


p/p SL and 
0.002376 sl/ff 


p = P a (psf)/RT (°K) 

d a 


i 










J 


M (ft/lb) 












Figure E-2. Flexibility Influence Coefficients, [*45/0] Test Specimen 





















Figure E-4. Flexibility Influence Coefficients, [+302/0] g Test Specimen 
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R » 3089.8 ft/sec 2 °K 

T (°K) * 0.5556 [ { ~ T - St °P ) + 460 ] 

a 2 

P a (psf) = Sg Hg u) H ^ 0 h Hg /12 = 70.72 h Hg (inches) 

where Sg„ = 13.6 
tig 

U true = 59 • 3 /h alc (inches)/a 


E.3 Experimental Tip Deflection Calibrations 


Calibration 1 raw data: 

tip lateral deflections (6 A ): ± 0.25 in; ± 0.5; ± 0.75; 

± 1.0 

obtain strain gauge readings from strip chart in units 
of mm 

bending strain gauge (channel 2) 
torsion strain gauge (channel 1) (0 c ) 


Calibration 2 raw data: 

tip angular deflections (e A >: ± 4.78 deg; ± 9.55 
obtain strain gauge readings from strip chart in units 
of nun 

By plotting each calibration point, obtain the following 


slopes: 
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a ll 

a 21 

a 12 

a 22 


5 c 

(cal 

1) 

5 A 

(cal 

1) 

9 c 

(cal 

1) 

5 A 

(cal 

1) 

5 c 

(cal 

2) 

0 A 

(cal 

2) 

9 c 

(cal 

2) 

9 A 

(cal 

2) 


From these, we can obtain the following transformation equa¬ 
tion: 


S(mm) 
c 

^ _ 

" a ll 

a 12 

I 6 A (in) 

0 (mm) 
c 


a 21 

a 22 

j 0 A (deg) 


Since we want the inverse relationship, the A-matrix is 
inverted to obtain the final expression. 



« [ A ] 


-1 









From this expression, we can get approximate values for the 
test specimen tip deflections knowing the strip chart strain 
gauge readings. 


I 


t 



EXPERIMENTAL NATURAL FREQUENCIES 
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LIST OF SYMBOLS 

A Amplitude of oscillation 

ag Two dimensional lift curve slope = 2ir 

b Semichord = c/2 

C(k) Theodorsen function 
Lift coefficient 
c Chord 

Flexibility influence coefficients 

Flexural modulus for an anisotropic plate 

E t Longitudinal modulus 

Li 

E m Transverse modulus 

T 

e Distance between elastic axis and aerodynamic center 

F fc Concentrated test load 

G t _ Shear modulus 

I4X 

g Structural damping coefficient 

i /=T 

k Reduced frequency = wb/U 

k.. Stiffness influence coefficient 
iD 

L_ Lift force at elastic axis 

E 

l Length of plate 

Aerodynamic Moment 
M fc Concentrated test moment 

m Mass per area 












Distributed lateral load 


Generalized force 
In-plane, on-axis lamina modulus 
Dynamic pressure 
Generalized displacement 

Kinetic energy 
Time 

Plate thickness 
Velocity 

Divergence velocity 
Flutter velocity 

Strain energy 

External work 

Lateral deflection 

Lateral deflection of elastic axis 

Tip lateral deflection 
P 

Complex eigenvalue = (1 + ig)/w 2 

Angle of attack 
Initial angle of attack 
Elastic twist 
, Tip angle of attack 

Mode shape (two dimensional) 

Ply angle 


Coupling factor 
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1 

» 


v__ Major Poisson's ratio 
v TL Minor Poisson's ratio 
Air density 

4 >^ Mode shape (one dimensional) 

a) Frequency of oscillation 
u) Torsional natural frequency 








